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Abstract 

Under the diffusion scaling and a scaling assumption on the microscopic component, a 
non-classical fluid dynamic system was derived in [5] that is related to the system of ghost 
effect derived in [4T] in different settings. This paper aims to justify this limit system for a 
non-trivial background profile with slab symmetry. The result reveals not only the diffusion 
phenomena in the temperature and density, but also the flow of higher order in Knudsen 
number due to the gradient of the temperature. Precisely, we show that the solution to the 
Boltzmann equation converges to a diffusion wave with decay rates in both Knudsen number 
and time. 
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1 Introduction 


Consider the Boltzmann equation with slab symmetry under the diffusive scaling 

ed t f + = ]Q(f, / £ ), (t, *,0 G K+ x R x M 3 . (1.1) 

Here f £ (t, x, £) > 0 is the distribution density of particles at (t, x) with velocity £, Q(f, /) is the 
collision operator which is a non-local bilinear operator in the velocity variable with a kernel 
determined by the physics of particle interaction. For monatomic gas, the rotational invariance 
of the particle leads to the collision operator Q(f, f ) as a bilinear collision operator in the form 
of, cf. [S] : 

Q(f, g m = ^ ( f (f(O9(£) + f(£)g(O-f(Z)9(Z*)-f(Z*)9(0)B(\Z-Um*dn, 

2 J R3 J §2_ V J 

with 0 being the angle between the relative velocity and the unit vector H. Here = {!! 6 
S 2 : (£ — £*) • H > 0}. The conservation of momentum and energy gives the following relation 
between velocities before and after collision: 

f£ = £-[(£-£*)■«] n, 

ui = & + [(£-eo-«] 

In this paper, we will consider the two basic models, i.e., the hard sphere model and the 
hard potential with angular cut-off, for which the collision kernel H(|£ — £*|, 6) takes the form of 


|£-&l, *) = !(£-£*,«)!> 


and 


= b(6)€L 1 ([0,n]), p>5, 


respectively. Here, p is the index in potential of the inverse power law that is proportional to 
r l ~P with r being the distance between two particles. 

Motivated by m, the following macroscopic and microscopic decomposition with scalings 
was introduced in [3] : 

f £ = ^d[p s ,£u e ,e s ] + £ G s ■ (1-2) 

Here M[ p s, eu E ,0 E ] is the local Maxwellian and G £ is the microscopic component. Moreover, the 
local Maxwellian Mye^ u efle\ is defined by the five conserved quantities, that is, the mass density 
p £ (t, x), momentum density m £ (t , x) = ep £ {t , x)u £ (t , x) and energy density e £ (f, x) + \\eu £ (t, x)\ 2 
given by 

p £ (t,x)= f f £ (t,x,£,)d£, 

Jr 3 

\{t,x)= [ i/}i(0f e (^x,0d4 for i = 1,2,3, 

Jr 3 

(t,x)= [ ip 4 ( 0 f e (t,x,£)d£, 


m- 


(1.3) 


P [e + yP 


as 


P £ (t,x ) 


M = M\ a e 1 (t, x, f) = — . 

[P ,£U Jv ' \J(2iTR0 £ (t, x)) 3 

Here '(p a (£ are the collision invariants: 


exp 


\^-£U £ (t,x)\^ 

2 R6 £ (t, x) 


(1.4) 


" V’oCO = i, 

< = & for i = 1,2,3, 

xo^iei 2 , 
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satisfying 


[ i’j{£)Q{h,g)d£ = 0, 
J R 3 


for j 


0 , 1 , 2 , 3 , 4 . 


Here, 9 £ is the temperature related to the internal energy e £ by e £ = § R6 £ with R being the 
gas constant, and eu £ is the bulk velocity. Note that even though u £ is of higher order, it is the 
scaled velocity that appears in the equations for the macroscopic variables p £ and 9 £ . 

The Boltzmann equation is a fundamental equation in statistical physics for rarefied gas 
which describes the time evolution of particle distribution. There has been tremendous progress 
on the mathematical theories for the Boltzmann equation with e being a fixed constant, such as 
the global existence of weak (renormalized) solution for large data in |11| and classical solutions 
as small perturbations of equilibrium states (Maxwellian) in |20l f3T %2] and the references 
therein, etc. 

On the other hand, the study on the hydrodynamic limit of Boltzmann equation is important 
and challenging. For this, it is well known that the classical works of Hilbert, Chapman-Enskog 
reveal the relation of the Boltzmann equation to the classical systems of fluid dynamics through 
asymptotic expansions with respect to the Knudsen number. For the hydrodynamic limit of 
Boltzmann equation to the compressible Euler system, we refer mm for the formal derivation. 
If the Euler system is assumed to have smooth solution, this hydrodynamic limit is proved 
rigorously in mm with and without initial layer respectively. 

However, it is well known that the compressible Euler system develops singularity in finite 
time even for sufficiently smooth initial data. The Riemann problem is the basic problem to 
the compressible Euler system, and its solution turns out to be fundamental in the theory 
of hyperbolic conservation laws because it not only captures the local and global behavior of 
solutions but also reveals the effect of nonlinearity in the structure of the solutions. There 
are three basic wave patterns for the Euler system, that is, shock wave, rarefaction wave, and 
contact discontinuity. For the hydrodynamic limit of the Boltzmann equation in the setting of 
Riemann solutions, we refer mm Eaiza m esj. 

Under the diffusive scaling, usually, the density function / e (t, x,£) is set as a perturbation 
of a global Maxwellian Mnou, i.e. 


f £ (t, x, f) = M [li0jl] + M [w] x, £) + • • • + £ n fn(t, x, £)) • 


(1.5) 


There has been extensive study on the hydrodynamic limit e —> 0 of the Boltzmann equation 
to the incompressible Navier-Stokes-Fourier system, for example, to justify the DiPerna-Lions’ 
renormalized solution in m of the Boltzmann equation to the Leray-Hopf weak solutions of the 
incompressible Navier-Stokes-Fourier system. For this, Bardos-Golse-Levermore [2] first studied 
this problem under certain a priori assumption. Recently, a breakthrough was achieved by Golse- 
Raymond in m which established a proof of such limit for certain class of collision kernels. After 
that, some progress was made for more general collision kernels, cf. |3fi| . In fact, there are also 
a lot of important contributions on this problem over the years, see USUI EB EU m S3 io] 
and the references therein. 

In the framework of classical solutions to the incompressible Navier-Stokes-Fourier system, 
it was proved in |9] that one can find a Boltzmann solution f £ (t,x ,£) such that /f is of order 
e 2 , but it is not clear about the amplitude /f at the initial time. Later, the Navier-Stokes- 
Fourier limit was proved for f £ ( 0,x,£) with small data in [3j. Recently, Guo in [2T] justified 


the diffusive expansion (1.5) when /i(0,x,^) has small amplitude while /f(0, x,£) can have 
arbitrarily large amplitude for i > 2 in a torus. This work was later generalized to some other 
settings, cf. [33, 28j[29]. Moreover, based on the L 2 — L°° estimate, Esposito-Guo-Kim-Marra 
m proved the hydrodynamic limit of the rescaled Boltzmann equation to the incompressible 
Navier-Stokes-Fourier system in a bounded domain if the initial data is small. 
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Notice that all the results under the diffusive scaling mentioned above are either about large 
perturbation of vacuum or small perturbation of a global Maxwellian. A natural question to 
ask is how about the perturbation of a non-trivial profile. The purpose of this paper is to study 
this problem in the setting of (1.2). 


In fact, under the assumption (1.2), when e —y 0, formally we have 


f £ — M[p e ,eu s ,0 s ] + £ G £ (1-6) 

which shows that in the macroscopic level, only the unknown limit functions p, 6 survive because 
the macroscopic velocity is zero. However, as shown in |3j and will be recalled in the next section, 
the equations of p and 6 are actually closely related to the scaled velocity u. Indeed, this diffusive 
scaling induces diffusion phenomenon for both the temperature 0 and density p, and the non-zero 
gradient of temperature induces a non-trivial flow in the higher order along the same direction. 
In this paper, we will construct such diffusion wave and study the hydrodynamic limit of the 
rescaled Boltzmann equation to such a diffusion wave global in time. 

The rest of the paper will be organized as follows. The construction of the diffusion wave 
and the main theorem will be given in the next section. We will reformulate the problem and 
derive some a priori estimates in Section 3. Based on the a priori estimates, the main theorem 
will be proved in Section 4. 


Notations: Throughout this paper, the positive generic constants that are independent of e 
are denoted by c, C, Ci(i = 1, 2, 3, • • •). And we will use || • || to denote the standard ^(M; dz ) 
norm, and || ■ || H i (i = 1, 2, 3, • • •) to denote the standard Sobolev P ? (M; dz) norm with z = x or 
y. Sometimes, we also use 0(1) to denote a uniform bounded constant independent of e. 


2 Construction of Profile and the Main Result 


We will drop the superscript e in the case of no confusion for simple notation. The inner product 
of h, g in L|(M 3 ) with respect to a given Maxwellian M is defined by: 

{h,g)M = [ Y~M0g(0d^ 

J R3 M 


when the integral is well defined. If M is the local Maxwellian M, with respect to this inner 
product, the macroscopic space is spanned by the following five pairwise orthogonal functions 


Xo(£) = ^rp M ' 

xM) ^ k mp M for i=1 - 2 ' 3 ’ 

(Xi,Xj) = hj, h3 = 0,1,2, 3,4. 


Using these functions, we define the macroscopic projection Pq and microscopic projection Pi 
as follows: 

' 4 

P 0 h = ^2(h,Xj)Xj, 

3=0 

k P\h = h — Poh. 

The projections Pq and Pi are orthogonal: 


PoP) = Po, Pi Pi = Pi • P(i Pi = PiPo = 0. 
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A function h(£) is called microscopic or non-fluid if 

J = 0, J = 0,1,2,3,4. 

Under this decomposition, the solution f(t, x, £) of the Boltzmann equations satisfies 

Pof = M, Pi/ = sG, 

and the Boltzmann equation becomes 

{eM + £~G)t + £i(M + eG) x = 2Q(M, G ) + eQ(G, G ), 

which is equivalent to the following fluid-type system for the fluid components (see [33] and 
for details): 


Ep t + {£pui) x = 0 , 

£(Epui)t + ( £ 2 pu\ + p) x = -e j t%G x d£, 

£{Epui)t + (e 2 puiUi) x = -e J ii£iG x d£, i = 2,3, 
e[p(e + + [epui(e + + epm] x = -£ J ^£i\€\ 2 G x d£, 

or more precisely, 

'spt + {epu{) x = 0 , 

e(epni) t + ( £ 2 pu\ +p) x = ^e(p(0)eu 1x ) x - e j £?© x d£, 

£{epui)t + (E 2 puiUi) x = £(p(9)£u ix ) x -£ £i£i@ x d£, i = 2,3, 

I i2 ^ 

I | o PI I ^ 

£[p(e + ^~)] x + [epui(e + + epui] x = e(k(9)9 x ) x 

4 3 f 1 

+ -e(e 2 p(0)uiui x ) x + ^2 £(£?p(0)uiu ix ) x - e / -fi|£| 2 <d x d£, 

- i =2 ' 

together with an equation for the non-fluid component G: 

£ 2 Gt + Pi(£iM x ) + ePi(^iG x ) = LmG + eQ(G, G ), 

where 


and 


G = L- M \P l {^ l M x )) + Q, 


0 = L-J^Gt + £Pi^iG x ) - eQ(G, G)). 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


Here Lm is the linearized operator of the collision operator with respect to the local Maxwellian 
M: 

L M h = Q(M, h) + Q(h, M), 

and the null space N of Lm is spanned by the macroscopic variables: 

Xj, J = 0,1,2, 3,4. 

Furthermore, there exists a positive constant ao(p,u,0) > 0 such that for any function h(£) E 
N- 1 -, see H91. 


(h,L M h) < -a 0 (v(\€\)h, h), 
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where z/(|£|) is the collision frequency. For the hard sphere and the hard potential with angular 
cut-off, the collision frequency z/(|£|) has the following property 

o < vq < kici) < K 1 + leiy 3 , 


for some positive constants vq, c and 0 < /3 < 1. 

In the above presentation, we normalize the gas constant R to be | for simplicity so that 
e = | R6 = 9 and p = Rp6 = \p9. Notice also that the viscosity coefficient p(9) > 0 and 
the heat conductivity coefficient n(9) > 0 are smooth functions of the temperature 9. And the 
following relation holds between these two functions, EKE], 

«(«) = f (2.4) 

after taking R = It should be pointed out that ( |2.4| ) is crucially used in the following analysis. 
In fact, in our analysis, it is required that 

inf k(9 ) > - sup p(9) 
e 4 q 


for all 9 under consideration. By (2.4), it is known that the above condition holds provided that 
the variation of the temperature is suitably small. 


Now we are in a position to derive the limit equations for (p,u,9) in the diffusive limit (1.6) 
formally. As [3], we assume that 

p £ = const + 0(l)e 2 , (2.5) 

then, as e —> 0, (2.2) x , (2.2 ) 2 and (2.2) 4 yields formally that 


p = const, 

Pt + ( pu\) x = 0, 

(pui)t + ( pu\) x + P* = p{9)u lx ) x , 
(p9) t + (pui9 +pui) x = {n{9)9 x ) x , 


( 2 . 6 ) 


where P* is unknown function. The equation (2.6) reveals how the zero order function p,9 
depend on the scaled velocity even though the macroscopic velocity tends to zero. 


With slab symmetry, in the macroscopic level, it is more convenient to rewrite the system 
by using the Lagrangian coordinates as in the study of conservation laws. That is, consider the 
coordinate transformation: 


U f(z.i) \ 

P{y,s)dy - (pui)(y,s)ds,s), 
(o,o) 7 


which is still denoted as ( x , t ) without confusion. Denote that v 
in the Lagrangian coordinates become 


j, the system (IT) and ( 2.1 ) 


and 


V V £ 

ev t - eui x = 0, 

£ 2 u u +Px = s J £lG x d£, 

£ £,iCiP'xd(,, i — 2,3, 

e(e + + (epui)x = -£ j 2 G x d£, 


(2.7) 


( 2 . 8 ) 
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respectively. Moreover, (2.2) and (2.3) take the form 


ev t - euix = 0, 

e 2 u u +p x = \e 2 {^-u lx ) x - e [ £i©i x d£, 

3 v J 

£ Uit = E ( — Uix)x £■ J £l£j©l:r^C i = 2,3, 

e(e + 7T- 1 —) t + (£pui) x = e(——0a,)® + -£'(——«iui x ) 3 

^ 'U O V 

+ ^£ 3 (^l Ui Ui x ) x - £ [ ^l\£\ 2 @lxd£, 

1=2 V ■' 


(2.9) 


and 


with 


and 


£ 2 G t - ^-^-G x + Pi(—M x ) + eP 1 (—G x ) = L m G + eQ(G, G), 

V V V 


G = Lyf(Pi(—M x )) + ©!, 

V 


®1 - l m {f'Gt -^ 


G x + -P 1 (£ 1 G x )-eQ(G,G)). 
v J 


The limiting equation (2.6) becomes 


( 2 . 10 ) 


( 2 . 11 ) 


p = const, 

Vt - U lx = 0, 

4 p{6) 

U lt P x — Q ( U lx)xi 
3 v 

@t T pu\ x — ( — 0x) x . 

U 


( 2 . 12 ) 


2.1 Construction of profile 


We will construct a background solution to (2.12) in this subsection. Without loss of generality, 
set 


20 2 
P= 3^ = 3’ 


that is 


(2.13) 


v = 6. 


(2.14) 


Assume the boundary conditions at the far fields given by 

6 6 

lim (v,0)(x,t) = (v±,6±), and — = — = 1, with 

X— >±00 V-L V— 


7^4 


(2.15) 


Note that if 0_ = 0 + , then v = 9 = 1 ,ui = 0 is a trivial solution to (2.12), and the diffusive 
limit of the rescaled Boltzmann equation to the incompressible Navier-Stokes-Fourier system is 
well studied as mentioned in the introduction. 

Noting (2.14), the equation (2.12) 4 is rewritten as 


0t + q Ulx = r~9x)a 

3 u 


(2.16) 






















Substituting (2.12) 2 into (2.16) and noting (2.14), we have the following scalar nonlinear diffusion 
equation 


9 t = ( a(9)9 x ) x , a(0) = 


3 k(6) 


, with lim 9(x,t) = 9±. 

00 x —>=too 


(2.17) 


From |Ij and m, it is known that the nonlinear diffusion equation ( |2.17[ ) admits a self-similar 
solution 0(rj) with r/ = satisfying the boundary conditions 9(±oo,t) = 9±. Furthermore, 


9(jj) is a monotonic function. Let 5 = \9+ — 9- 1, then 9{t,x ) has the property that 

.2 _ 

Too. 


9 x {t,x) = t^=}Le 4 “i # ±)d+‘) 


as x 


Define 


\/1 +t 

v,ui,0) = (9,a(9)9 x ,0)(x,t) 


(2.18) 

(2.19) 


then it is easy to check that (v,ui,9) satisfying (2.12) as 

‘ _ 29 2 

P= 3r, = 3' 

Vt - Uix = 0 , 

4./i(0). . 

Hit + P x ~ n (-~— U lx)x, 

3 v 

9t "F puix — ( ~ 9x)xi 


( 2 . 20 ) 
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where P* = — a{9)9 t + —-p-{a(9)9 x ) 


Remark 2.1 By (2.19) and (2.20), we actually construct a diffusion wave to the limit system. 
On the other hand, if 9- < 9+, then ui = a(9)9 x > 0, that is, the variation of temperature 
along the x-axis induces a nontrivial scaled flow along the same direction, see Figure 1. The 
case 6L > 9 + is similar, see Figure 2. 




Remark 2.2 The construction of the profile (v,ui,9) is motivated by the viscous contact wave 
of compressible Navier-Stokes equations, see JUJ/, l2Uf and The viscous contact wave is 

used to approximate the contact discontinuity for compressible Euler equation and its pressure 
keeps constant. 


In order to justify the hydrodynamic limit of the rescaled Boltzmann equation to the limit 
system (2.20), if we use the profile (v,ui,0), then some non-integrable error terms with respect 
to time coming from the non-fluid component for the system about perturbation. Therefore, 
one needs to construct another profile (v, eu , 9) for the rescaled Boltzmann equation, based on 
(v,ui,9). For this, we require that the approximate pressure p satisfies 


V = | + 0( l)e 2 =p+ + 0(l)e 2 . 

6v 3 


( 2 . 21 ) 
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Motivating by [27j, we first notice that the main part of the non-fluid component in the solution 
G and part of 0i defined in (2.11), are given by 

w = ~ P M (-^l(£l-^x)) = J^ V Q P M {-^l[Cl( + C ■ £'U x )M]}, 


and 


©i = L m {~ p i{ii w x) - eQ(w,w)), 


respectively. To distinguish the leading term coming from the non-fluid component, we rewrite 
the Boltzmann equation (2.9) as 


' ev t - £u lx = 0, 

2 

£ 2 uit+p x = \e 2 (^-u lx ) x - [ ^©i x d£, 

V 2 J=1 J 

£ 2 u it = £ 2 {^-u ix ) x - f £iti@{ x d£, i = 2,3, 

V 3=1 J 

e(e + ^-) t + (£pui)x = £(^p~0 x ) x - j ^i|^| 2 ©L^ + H x, 

3=1 


( 2 . 22 ) 


with 


e^G t -L M G=-—P 1 M 


Rv9 


1C ~ eu\ 
29 


■ (6 - 0) x + £ • (eu - £u) x )M] 


£ 2 U\ 


+ - ~G X - -Pi(6Gx) + eQ(G, G) - £ 2 G u 

V V 


(2.23) 


where 


G = ikeLuiPiU^Ox + C • eu x )M]}, G = G — G, 




H = ^j-^P-UiU 1x + Yli=2 £S ^r U i U iX’ 

*1 _ r- 1 ( e 


3 .3 M0, 


Q\ = L^[^P 1 (CiG x )-£Q(G,G) , 


0f = L-^{fG t - ^G x + zPifaG x ) - £Q(G,G ) - 2 eQ(G,G)), 


satisfying 


J2G{ = Gi = LM 1 (e?G t 
3=1 


—G x + -Pi(£iG x ) - £Q(G, G)). 

V V 


(2.24) 


Here, the function (v,£u,0)(x,t) is the profile to be constructed. 

Since the velocity £u decays faster than (v, 6) in time, the leading terms in the energy 
equation (2.22 ) 4 are 


£0 t + £pu lx 




— £ 


^ilci 2 ©LrfC- 


(2.25) 


By the definition of 0}, it holds that 
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-e/^i|£| 2 ep£ = e 2 JVi + e 3 Fi, 

Nl = fll9 x 9 x + fl2V x 9 x + f\'.$x + /l 4 ® XX ; 

k 1-^1 I = 0(1)[(|'U 3 ;| “I - |^a;| “I - \@x\ T e|Wx-| "I" | ) |ltl I “I - |'Ux^o;| “I - I'U'Xxj]) 


(2.26) 


where the coefficients fij,j = 1, 2, 3,4 are smooth functions of (v, eu, 0). By (2.21), it is expected 
that the profile (v, eu , 9) for the Boltzmann equation satisfies 6 = v. Thus, by choosing only the 
leading term in (2.25), one obtains that 


eOt — e(a(9)9 x ) x H— —N± x , 
5 


(2.27) 


where a{6) is given in (2.17). Thus the leading part of (2.27) is the nonlinear diffusion equation 
(2.17) and an explicit solution 0(- 7 =) is given with the boundary conditions Q(±oo,t) = 9±. 

To include more microscopic effect, let the profile 0 & 9{ v /j^ ) + e0 n * (x,t), where 9 n *{x,t) 
represents the part of the nonlinear diffusion wave coming from the non-fluid component. More¬ 
over, the term 9 n f (. x,t ) in the form of - 7 =T>i(- 7 =) is from Af in (2.27). Note that 9 n f (. x,t ) 

decays faster than 9(x,t ) so that it can be viewed as a perturbation around profile 9{x,t). To 
construct 9 n f ( x , t), we linearize the equation (2.27) around 0(x, t ) and keep only the linear terms. 
This leads to a linear equation for 9 n f (x,t) from ( 2.27[) 


Of = {a{9)9 n J) x + (a\9)§ x 9 n % + -N lx , 
where Ni = (/n + /i 2 + j\:i)(9 x ) 2 + fu9 xx with fy = fy(v, 0, 9),j = 1, 2,3,4. Let 


(2.28) 


gi(x,t) = I 9 nf (x,t)dx, 

J — OO 


then integrating (2.28) with respect to x yields that 

„ „ „ 3 . 

git = a(9)gi xx + a! (6)0 x gi x + -Ni. 

5 

Note that N\ takes the form of Jl+t ) an< ^ sa tishes 


(2.29) 


llVil = 0(l)S(l + t)~ 1 e 


as x 


Too. 


We can check that there exists a self-similar solution gi(g),g = x/y /1 +1 for (2.29) with the 
boundary condition gi(—oo,t) = 0, gi(+oo,t) = 5\. Here <5i satisfies 0 < <5i < <5. Note that even 
though the function g\(x,t) depends on the constant 5i, 9 n f(x,t) = g\ x {x,t) —> 0 as x — > Too. 
That is, the choice of the constant <5i has no influence on the ansantz as long as |5i| <6. From 
now on, we fix 5\ so that the function g\(x, t ) is uniquely determined and its derivative g\ x = 9 n * 
has the property 

i - 2 

\9 nf \ = \gi x \ = 0(5){1 + t)~2e Me±)W , as i-> Too. 

Now we follow the same procedure to construct the second and third components of the velocity 
profile denoted by eut, i = 2,3. That is, the leading part of the equation for eu t coming from 
(pT22l) is 


e u it = e ( 


2/0(0) 


Ur 


— £ 




(2.30) 
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For i = 2,3, one gets 


'-£ f £i£i@\d£ = e 2 Ni + £ 3 Fi, 

F/ — fn9 x 9 x “I - fi 2'Vx@x T f a 0 x F fa9 XX i 

^ | Fj | — 0(1) ((I'U.j | + | O x | T \9 X \ F £|'Mx| F sI'W.t |) \v>x | “1“ |^a:||^a;| “1“ |’U , £a;|)i 


(2.31) 


with smooth functions fij,i = 2,3, j = 1,2,3,4. Notice that the symbols N t and F), i = 2,3, 
used here are for the convenience of notations. 

From (2.30) and (2.31), we expect that the profile Ui(x,t ) takes the form of ) 

and satisfies the following linear equation 


E 2 u it = £ 2 {^p-u ix ) x F £ 2 Nix , i = 2, 3, 

6 

where iV* = (fn F fi 2 + fi3){9 x ) 2 F fnOxx, fij = fij(v, 0, 0),i = 2,3, j = 1, 2, 3, 4. 
Denote 

rx 

gi(x,t) = / Ui(x,t)dx, 


(2.32) 


then integrating (2.32) with respect to x, one has 


9it — ~ fjixx F Ni. 

6 


(2.33) 


For given 9 , we can check that there exists a self-similar solution gi(g) with g = ^== with 
the boundary conditions gi(—oo,t) = 9 i(Foo,t) = 6i, where 5i satisfies 0 < <5* < 5. As we 
explained before, the choice of the constant <5j is not essential. From (2.18), we fix 5i so that the 
function gi(x,t ) is uniquely determined and the derivative gi x = iii (i = 2,3) has the following 
property 


\eui\ = \egix\ = 0(l)5e{l F t) 4( >( e ±x i + t ) ; as x -> ±oo, 


where b(9±) = max{a(0±), 

In summary, one can define the profile (u, eu, 6) for the Boltzmann equation as follows. To 
satisfy the conservation of mass, one needs 


ev t - euix = 0. 


By letting v = 9 + eO n f , one gets 


3e 2 


eu\ = £[a{9% F £a(0)9f F £a\9)9 x 9 nf ) + —N x . 

o 


(2.34) 


However, by plugging (2.34) into the momentum equation of (2.22), we have a non-conservative 
term containing £ 2 N\t■ To avoid this, one defines 

£U\ = £[a(9)9 x F £a(9)9™f F £o!(9)6 x 9 n ^]. 

Similarly, to avoid the non-conservative term (\u\ 2 )t in the energy equation, set 

0 = 0 ns + £0 nf - -\£U\ 2 . 

2 
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Therefore, the profile (v,eu,6) is finally defined as: 
v = 9 + eO nf , 

eu\ = e[a(9)6 x + £a{9)9 x * + ea'(9)9 x 9 n ^], 
£Ui = £(Jixi i = 2, 3, 

9 = 9 + £9 n f — ^|efi| 2 , 


(2.35) 


where 9 is given by (2.17), 9 n f by (2.28) and g ,, % = 2,3 by (2.33). Then a direct but tedious 
computation shows that 

3e 2 - 

£V t - £Ui x = —N lx , 

5 

2 - - 4e 2 ,n(9) _ - 

£ 'U'\t Px — Q ( Z 'U'lxjx H - R\xi 

6 v 


e 2 Uit — £ 2 (^—+^-Uix) x + £ 2 Ni x + Ri x , i — 2, 3, 

I -12 V 
\£U' 


e(e + 1 i ) f + {£pu 1 ) 3 ; — e( ^ 9 X ) X + H x + s 2 N\ x -— N\ x + R^x, 


(2.36) 


where 


Ri = e 2 [a(9)9 t + (, a(9)9 n f ) t ] +p-p + - ^(^W) 

X 2 o V 

= 0 (l)cfe 2 (l + t) _1 e 4c ( fl ±)( 1+t ), as x — > ± 00 , 
g(9) n(S) 


Ri = £ 


-]£Uix + £■ (Ni — Ni ) 


' 9 v x2 

= 0(l)<5e 3 (l + t)- 3 / 2 e _ 4=(«±)(i+‘), as x —> ± 00 , i = 2,3, 

R, = \le{a{9)9 x + a(9)9 n J + a’(9)9 x 9 n f) - e^9 x 

L3 A _ _ v 

+{p - P+)eu\ + £ 2 {Ni - Ni) - H 

s 2 

= 0(5)e 3 (1 + t)~ 3 / 2 e 4c ( s ±)( 1+t ), as x — > ± 00 , 


(2.37) 


(2.38) 


(2.39) 


Ni = 0(1)5(1 + t) l e 4a («±)( 1 + t ) ; as x -» ± 00 , i = 1,2,3, 


(2.40) 


with c(9±) = max{a(0±), \b{9±)}, Ni, i = 1, 2, 3, and H are the corresponding functions defined 
in ( |2.24 ), (2.26) and (2.31) by substituting the variable (v,£u,9) by the profile (v,eu,9). Note 
that the decay rate of Ri, i = 2, 3,4 is of order e 3 (l+t)~ 3 / 2 . Furthermore, even though the decay 
rate of R\ is still e 2 (l + t) _1 , it is sufficient to obtain the desired a priori estimates through some 
subtle analysis coming from the intrinsic dissipation mechanism in the momentum equations as 
shown in the following. 

Define 

|£ - £u\ 2 \ - _1 (\ ^ 


M = 


v 1 


y / ( 2irR9 ) 3 


exp 


and 


2R9 
f = M + £G 0 . 


. G„ = L-J ( ) , 


Then it follows from (2.36) that 


sft - —fx + -Cifx — L m Gq + eQ{Gq, Go) + Rf, 

v 


v 


(2.41) 
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where 


Rf = t, i)M + e 2 G ot - e £ -^G 0x + eP\ (^G 0x ) - eQ(G 0 , G 0 ), 

J v \v / 


and \B 2 (x, t, £)| = 0(1)<5(1 + t) ae 4c (®±)( 1 + t ) |£| 3 , as x -» ±oo. 


Remark 2.3 From the definition of [v,u\,9) in (2.19) and the definition of (v,u\,9) in (2.35), 
it holds that 


\(v — v,ui — iii,0 — 0)(x,t)\ = 0(l)5e(l +1) a e 4c («±)( 1 + t ) ; 
that implies that the ansantz (y,ui,9) well approximates (v,u\,9) when e is small. 

2.2 Main result 


(2.42) 


Now we consider the system rt2.9|-(2.10) with the initial data 


(v, u, 9) |t =0 = (v, u, 0)(x, 0), G(x, t)\ t =o = G(x, 0). 

Then the main result in this paper can be stated as follows. 


(2.43) 


Theorem 2.4 Let ( v,u,9)(x,t ) be the profile defined in (2.35) with strength 5 = |0+ — |. 


Then there exist small positive constants 5 q and £o and a global Maxwellian M* = 

such that when 5 < 6 q and e < £q, the Cauchy problem (2.9)- (2.10) with the initial data (2.43) 


has a unique global solution ( v,u,9,G ) satisfying, for any sufficiently small but fixed positive 
constant 9 > 0 , 

(v — v, eu — £u, 9 — 9)(t )||| 2 < C'v / 5e 3 (l + t)~ 1+c °^ , 

(v-V,£U-£U,0-9) x {t) f Ll < CV5£ 2 (l+t)-l + » +c °' r5 , 

fxx{t)\\ i 2 (i 2 (_i_)) + \\(v-v,eu-eu,9-9) xx (t)\\ 2 L 2 < C\f5(l + ^ ^ 

(G - G){t)\\ 2 L ,2l(^)) ^ C ^ 1 + 


(G G) xW || i2 (L 2 (7 L }) 


< CVS( 1 + t)“i+ ,? + c °' / *, 


that implies that 


I || (v — V,£U — £U, 9 — 9) (t)||£,oo < C5i£ 4 (1 + t) 8 + 4^, 
[\\(v - V,£U- £U,0 - 9) x (t)\\ L ™ < C6*£5(l + t) _ 3 +1? , 

where C and Cq are positive constants independent of £ and 5. 


(2.45) 


The following result justifies the hydrodynamic limit of the rescaled Boltzmann equation 


(1.1) to the diffusion wave (v, u\, 9) global in time. 


Corollary 2.5 Under the conditions of Theorem \2.J\ from (2.42) and (2.45), it holds that 
f\(v-v,9-0)(x,t)\<Ce(l + t)-z->0, 

< i i as £ —> 0, 2.4b 

1 \(ui - ui)(x,t)\ < C£4(l + t) 2,-^0, 


that is, the fluid part (v, u\, 9) of the solution of the rescaled Boltzmann equation (1.1) converges 


to the diffusion wave solution {y,u\,9) of (2.20) in the sense of (2.46) as £ —> 0, which reveals 
that v and 9 are diffusive. 
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Remark 2.6 The above Corollary shows that if the zero order function in (1.5) is not a global 
Maxwellian, then one has to consider the effect of diffusive wave in the diffusive limit of rescaled 
Boltzmann equation (1.1). 

Since the scaled velocity u\ is actually induced by the variation of temperature 9 , i.e., 
u\ = a(9)9 x . The following result shows that the scaled velocity u\ is also induced by the 
variation of temperature 9 in some sense when e is small. From the definition of 9{rf) with 
rj = in (2.17) and (2.18), it can be seen that 9 is monotonic. To be definite and without 

loss of generality, let us assume that that is, 9 is monotonically increasing. Then there 

exists a positive constant r/o > 0 such that 


9'{rf) > crj 0 5 , for \rj\ < rj 0 , 


(2.47) 


where c Vo depends on r/o and c Vo —> 0 as rjo —> +oo 


Corollary 2.7 Under the conditions of Theorem 2-4 and 9- < 9 + , for any fixed r/o > 0, there 
exists a small positive constant E\ = £1(770) < £0; such that if e < e\, then it follows from (2.47) 
and (2.45) that 


0 < CiVi+t K ~ 

0 < \9 X < 9 x (x, t) < M x , 


for |x| < 770(1 + t) 2 , t > 0, 


(2.48) 


that is 


3Ci 


6 x (x,t) < ui(x,t) < 2Ci6 x (x,t), for |x| < ? 7 o(l + t) 2 , t> 0, 


(2.49) 


where C\ is a suitably large positive constant depending only on 9±. In particular, (2.49) implies 
that variation of the temperature induces a non-trivial flow of higher order in the following 
parabolic region 

j(x,f) : \x\ < 770(1 + t) 2 , t > 0 |. 


3 Stability Analysis 


In this section, we will investigate the stability of the profile constructed in (2.36) for the 


Boltzmann equation (1.1). This section is organized as follows: in Section 3.1, the fluid type 


system (2.2) is reformulated in terms of the integrated variables; Section 3.2 is devoted to the 


lower order estimate, while Section 3.3 is for the derivative estimate. 

3.1 Reformulated system 

We now reformulate the system by introducing a scaling for the independent variables. Set 


x t 

y = T = ~o- 

tr 


(3.1) 


In the following, we will also use the notations (v, u, 9 )(t, y) and (v, u, 9 )(t, y ), etc., in the scaled 
independent variables. Set the perturbation around the profile ( v,u,9)(r,y) by 


= v — v, = eu — eu, Q = 9 — 9, 
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and 


(M,iy)(y,T) = 


f 


\£U\ 


£U\ 


> V’j H ) — (^ + 2 


(z, r)dz. 


Then we have ( 0, -0) = (4>, T) y and ( + || T y | 2 + X^=i £u i^iy = W'r 

Subtracting (2.36[) from the equation (2.22) and integrating the reduced system yield 


2 n 

- Tiy = - - 2 


5 P+ 


£ Ni, 


^lr + 


4e //i(0) 4 /i(0) _ \ f y2r\j ta 

P“P = ~~ 3 ~ uiy ) _£ S J ^i 0 i^ 

i=i 


— i?i, 




— L u iy - -zr Lu iyj+£ 2, {N i -Ni)+£ i F i -£ j 6 ?i©i^ - Ri,i = 2,3, 

IU T + eptti — ephi = 


(3.2) 


V) + (ff - ff) + e 2 (iVi - IVi) + £ 3 Ti 


w y ~ w y 

V V 


-e Jl^el 


2 di - R 4 + -£ 2 1Vi. 
5 


Since the variable IF is the anti-derivative of the total energy, not the temperature, it is more 
convenient to introduce another variable 


It follows that 


W = W-eih9 1 . 


( = W y - Y, with Y = -|F y | 2 - + £u 2 ^ 2y + £u 3 ^ 3y . 


Using the new variable W and linearizing the left hand side of the system (3.2) by using the 
formula of H in (2.24) give that 


^r-^iy = -|e 2 Wi, 

,Tr P+ * , 2 m , 4 ^( 0 ) 

^lr _ “1“ Q _Wy Q _ ^lyy o ( - J^^ly 

v 3v 3 v 3 v v 


— £ 


J2 [+ Ji + l=Y -Ri = + Gi, 


3=1 


3 v 


ViT = —4' iyy + e(— - ^-)u iy + £ 2 (Ni - N t ) + £ 3 F; 


— £ 


J 


V 

2 — Ri = ^ 4 /%yy + Qii i = 2,3, 

v 


(3.3) 


Mr + P+* i„ = yy W ro + (yy - ^p-)0y + e 2 (^i - JVl) + + y i„ 

3 — 

+ _ £ ^1 UiU iy ] - £Uir^i + J 2 - £ f ^£l|£| 2 ©l$; 


1=2 


+ e 2 «i Z)Li I C 2 @i^ - + |e 2 ^i + ewi-Ri - 


J 3 


= ^-Wyy+l^m + Ql, 
v o 
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where 


Ji = P ~^ ± %-\p-P + ^y-^ ; (0-e)\ = 0( ml + {o- of + Ml 4 ), 
J 2 = (p+ -p)^ly = 0(l)^l + 4 >l y + ( 0 - 6) 2 + Ml 4 ), 

Qi = - nr)«i» - * E [am + J, + |r - a,, 


3 v v 


Qi = e( 


/M) /M) 


1=1 


V V 

,k{9) k(0 ). 


)u iy + e 2 (Ni - Ni ) + e 3 Fj - e / — -R,, '« = 2 ,3, 


(3.4) 


.2/ AT AT \ | -.3 771 , 4 £ //(0) 


Q4 — (- —)®y + e Ml ~ M) + £ + 

v v 3 v 

+M'— 


ly 


i=2 


H iy 

V V 


UiUiy\ - eu lT ^i + J 2 - £ J ^£i|£| 2 0?d£ 


^ — 

+ £ 2 h! ^ [ gm - ^Yy + e«i^i - M 


i=i 


The equation of microscopic component G given in (2.23) in the coordinate (y,r) becomes 

~ ~ 1 If — £771 1 1 

vG t - vL M G = -—P^lL—L-Sy + z. -tpy'jM] 


+ £U\Gy — vP\(fi\Gy) + £vQ(G, G) — vG r 


In the scaling of (3.1), the equation (2.7) reads 


fr ~~ e ~fy + ~fy ~ zLmG + £ 2 Q(G, G). 
v v 


Set 


/ = /-/, 


(3.5) 

(3.6) 

(3.7) 


then from (3.6) and (2.41), we have 


vf T — suify + £if y — £VL M G + £[vLmG — vL m Gq] + £ 2 \yQ{G, G ) — vQ(Gq, Go)] 

- Mr + iffy - evRj. (3.8) 

Note that to prove the main theorem in this paper, it is sufficient to prove the following a 
priori estimate in the scaled independent variables based on the construction of the approximate 
profile. 

Theorem 3.1 (A priori estimate) For any sufficiently small and fixed positive constant id > 
0, there exist small positive constants 6 2 > 0,£2 > 0 and a global Maxwellian M * = 
such that if 5 < 5 2 and £ < £ 2 , then the Cauchy problem (3.3), (3.5) and (3.8) admits a unique 
smooth solution satisfying 

' IIM , W0(r)||£oc < CM, IIM if, C)(r)III, < cM 2 (i + £ 2 t)- 1+Co ^ 5 , 

^,C)y(r )\\ L 2 + IMMCMMlUg + Y L L ^JI^ d ^dy - CVS£ 3 (1 +£ 2 r)-i +1?+Co ^, 

lan 12 


l «[=2 


IM 

M* 


d^dy < C\f5{\ + £ 2 t) 2 , y] 


l«l=l 


YY-d^dy < CVS£( 1 + £ 2 t )-1+^+ c oV5^ 

lVl-jf 

(3.9) 


where C,Cq are positive constants independent of 6 and £. 
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In the next subsection, we will work on the reformulated system (3.3) and (3.8). Since the 
local existence of the solution can be proved similarly as the discussion in m and [42] , we will 
omit it here for brevity. To prove the global existence, it is sufficient to close the following a 
priori estimate: 


N(t)= sup [e 1 ||(<h,T,ir)||ioo +e 2 ||(0,V',C)ll|2+e 3 ||(</V ifiy, Cy)\\h 

(3 ' 10) 


0<S<T 

+ 1 


AT 


|o|=l 


M=2 


where Ao is positive small constant depending on the initial data and AT is a global Maxwellian 
to be chosen later. 

Before proving the a priori estimate (3.10), we list some lemmas based on the celebrated 
H-theorem for later use. The first one is from [El- 

Lemma 3.2 There exists a positive constant C > 0 such that 

r vmr l QUi9? 


M 


-d£ < C 


^ X 5 


£ + [ fix 

M ^ L 3 m L 3 m ^ 


M 


}> 


where M can he any Maxwellian so that the above integrals are well defined. 


Based on Lemma 3.2 the following three lemmas are from m- 


Lemma 3.3 If 6/2 < 6* < 6, then there exist two positive constants a = a(p, u , 0; p*, it*, Of) > 0 
and rjo = rjo(p, u, 6; p*,u*, Of) > 0 such that if \p — p*\ + \su — uf + \6 — Of < rjo, we have for 


OLmO 

AT 


> a 


AT 


dfi, 


where AT = and the definition of M = Mu £U m can be found in (1.4). 


Lemma 3.4 Under the assumptions in Lemma 3.3, we have 

vmr l h 2 


hjf \L-Jh\ 2 d( < a- 2 

m) \L-Jh\ 2 d(<i >- 2 


M 

mr'h 2 


df, 

df 


. JR 3 AT Jk3 AT 

for each h(£) G N 3 -. 

Lemma 3.5 Under the conditions in Lemma 1 there exists a constant C > 0 such that for 
positive constants k and X, we have 

, f giPim k 92) 


AT 


-df- 


, c r + 


AT 


AT 


Note that (3.10) also gives the a priori estimates on \\((f> T , ip T , C T )I|, ||d a (<TV’>C)ll and 


\d a G\ 2 
~M~, 


-dfdx 


(|a| = 2). In fact, from (3.1), (2.8) and (3.10), one has 
||(T-,Vv,Ct)|| 2 < C\\(v T ,£u T ,6 T )\\ 2 + C5£ 3 (l + e 2 r)"i 

< c[\\{p y -py, spuiy - £puiy)\\ 2 + \\(py,£pu ly )\\ 2 + £ 2 J J ^ ^ d^dy^j + Cd£ 3 ( 1 + e 2 r)“i 


< C/Wi^y^yXy^l 2 + d£ 2 U(f,fi,C )\\ 2 + £ 2 

< C(<5 + A 5 )e 3 , 


I Gy 


AT 


» 

-d£dy) +Cd£ 3 {l +£ 2 t)~I 


(3.11) 
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where we have used the fact that 


f ( J tlGydtfdy <cj f ^ didy. 


To derive the a priori assumption on 


, ij), Oil, (|a| = 2), we use the definition of p, m = epu 


and p{6 + ^\eu\ 2 ). Let |a| = 2, by (1.3), one can obtain 


Y (p,m,p(0 + -M 2 ))|| 2 < C 


\d a f \ 2 
AL 


-didy 


< C 


-d£ + C5e 3 ( 1 + t )~5 < C(\ 2 0 + 6)e 3 . 


AL 


(3.12) 


This yields that 


E W^AX)\\ 2 <C [ -J^ d ^dy+ E \\d P (^X)\\ 2 + C5e 3 (l+t)-i <C{\l + 5)e 3 . 

I« 2 J J * |/3|=1 

(3T3) 

Finally, one has 


|d“G | 2 


< C 


M, 
\d a f\ 2 
AL 


d£,dy < C 


\d a f\ 2 


AT 


-d£dy + C 


\d a (M — M)\‘ 


AL 


d£dy 


(3-14) 


d&y + C E 

M=l,2 


,^,Olr+ C'fe 4 (l + t)- 2 < C(eo + 5)V, |a| =2. 


3.2 Lower order estimate 


We are now ready to derive the lower order estimate. Multiplying (3.3) 1 by p+4>, (3.3) 2 by hihi, 
(3.3 ) 3 by T,;, (3.3 ) 4 by g^-VF with p+ = | respectively and adding all the equations, one can 
obtain 


+ tJ—W 2 + ^T 2 + J E 


3p+ 


+ + e + 2 ^w 2 


i =2 


i=2 


2>p + v 


y 


= ^e 2 JVi(—$ + ^-W) + hyT 2 + vQi®i + V ^ WQ 4 
5 3p+ 2 “ 3p+ 


(^)y*i*u, - E ~ + (■ • • V 

i=2 ^ + 


(3.15) 


Here and in the sequel the notation (• • • ) y represents the term in the conservative form so that 
it vanishes after integration. Since it has no effect on the energy estimates, we do not write 
them out in detail. 

Note that the term Qi^i contains (l + t) -1 \I , i which can not be controlled by the dissipation 
from the viscosity and heat conductivity. So is the term fVi(—<h + ^~W). As we will see 
later, an intrinsic dissipation associated with the profile is derived by the diagonal method and 
weighted energy estimate to control the above two terms. Let us consider the equations for the 
conservation of the mass, the first component of velocity and energy by defining 


v = ($,*!, wo*, 


where (-,-,•)* means the transpose of the vector Then from (3.3), we have 


V T + A\Vy - A2Vyy + A 3 


(3.16) 
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where 


/ 0 

-1 

0 \ 


( 

0 

0 

0 \ 

1 _p± 

1 v 

0 

§ 5 ) 

A 2 = 


0 

M9) 

3v 

0 

\ 0 

p+ 

0 / 



0 

0 

m l 

v / 


^3 — ( _ -£ 2 -Ni, Qii Qa + -£ 2 NiY. 

0 o 

Direct computation shows that the eigenvalues of the matrix A\ are Ai, 0, A 3 . Here A 3 = — Ai = 
A. The corresponding normalized left and right eigenvectors can be chosen as 

h = h = >^75(1,0,—), h = V^7l0(-1, 

3A 3 3 p + p + 3A 3 3p+ 


n = V WiO(-l,-A 3 ,p+) t , r 2 = V2/5(l,0,-p + ) t , r 3 = y/3/W(-l, A 3 , p+f, 


such that 


with 

Let 


liTj — 5ij , z,j — 1,2,3, LA\R — A — 

L = (h, l 2 , h)\ R = (ri,r 2 , r 3 ). 
B = LV = (h,b 2,63), 



then multiplying the equations (3.16) by the matrix L yields that 

B r + A By = LA 2 RByy + 2LA 2 RyBy + [{L T + A Ly)R + LA 2 Ryy\B + LA3. 


(3.17) 


A direct computation shows that LA 2 R = A4 is a non-negative matrix. From (3.17), we will 
apply weighted energy method to derive an intrinsic dissipation. Since we have assumed that 
9 y > 0. Let v\ = then |ui — 1| < CS. Multiplying (3.17) by B = (u^hi, 62 , uf ” 63 ) with a 
large positive integer n which will be chosen later, we have 


\v^b\ + \b\ + K n& i) r - Cir)r b 1 - ( V)- 6 i + ByA^By + BA 4y B y 
-^^(nXiVxy + vx\i y )b\ + \v^-\n\ z v ly - uiA 3y )bl 
= 2BLA2R y B y + B\LiR H- LA2Ryy\B + BAL X RB + BLA% + (• • • )#. 


(3.18) 


Let 


Ei = 

Ki = 
Note that 


3,2 + Ye' 2+ + I'E ' t ') dv + / ( f 65 + \ bl + V \ ibl)dv ' 


P+ a .2 
2 


i=2 


' 3 ly 


'E—K + 2 ^W] + By A 1 B y )i«. 


i=2 


3 p+v 


(B — B) y A±B y dy < C5 J \B y \ 2 dy + C6~ l J \8 y \ 2 \B\ 2 dy 

<C£ 2 5(l + t)~ 1 E 1 +C5Ki + C5 J \$ y \ 2 dy. 


(3.19) 
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Similarly, the terms in the last second line of (3.18), BA^ y B y , BLA 2 R y B y and B[L T R + 
LA 2 R yy \B satisfy the same estimate. For BAL y RB and BLA 3 , we need to use the explicit 
presentation. By the choice of the characteristic matrix L and R, we have 


A L V R = 



Thus 




LA, = 


^(^i+q 4 )-v§s 


_ r^Qx \ 


2 Q 4 
5 p+ 


\ V 1 a ^( £2N 1 + Q4 ) + \Jl% ) 


where 


BAL y RB = -A 3 y (v?bi + - vfhh - v^ n b z 3 ), 

BLA 3 = J £ 2 Ni(v™bi + v^ n b 3 ) + qiv%bi + q 2 b 2 + q 3 v\ n b 3 , 
V 15 p+ 


/ 2 1 5 Qi / 2 Q 4 /2 1 /5Qi 

91 = V ^ Q4 " V 6 V 92 = V 5^’ 93 = Vl^ 4+ V 6 V 


Combine (3.15), (3.18), (3.19)-(3.20), we have by choosing n sufficiently large, 


1 


(3.20) 


£ 1t + -Ki + 2 / + ^)|dy < Ce^(l + t)~\E 1 + 1 ) + CS $ 2 y dy + I nf , (3.21) 


where 
Rf = 


-)Qi^>idy+ / Qj'S'jdy + 


i =2 


3p+ 


WQ±dy + / (giui&i + ? 2&2 + n b 3 )dy.{3.22) 


Here we have used the fact that 


- <h + -—W = y/5/6(bi + b 3 ), 

3 P+ 


(3.23) 


and 


e 2 J |iV 1 |(| 6 1 | + \b 3 \)dy < C6 j \9 y \(b 2 + b 2 3 )\dy + Ce 2 5(l + t)~\ 
and for n large enough, 

~\ v i ^{nXiUiy + 2v\Xi y )b\ + ^v 3 n ~ l (nX 3 viy - 2vi\ 3y )b 2 3 - BAL y RB > 3\9 y \(bl + &§)• 


Even though Q± contains the term R.\ with the decay rate , the terms in (3.22) involving 
Q 1 have factor b\ or b 3 because 


Ti = \/ 3/IOA 3 (63 — 61 ). 


(3.24) 


Thus the terms vQi^i, q\v™b\ and q 3 v 1 n b 3 can be controlled by the intrinsic dissipation on 
b\ and b 3 as shown later. The estimates on the other terms involving Qi (i = 2,3,4) are 


straightforward because from (2.37)-(2.40) and (3.4), they decay at least in the order of e 3 (l + 


t) 3 / 2 . For brevity, we only estimate / vQi^idy and / f/ 2 ^ 2 ^ 2 / as follows for illustration. 

Estimate on / vQ\^\dy: 



















21 


From (3.24), we have 


vQi'&idy = 


^ / vQi\ 3 (b 3 - b\)dy. 


(3.25) 


Here we only consider the integral 


h= vQi\ 3 bidy , 


and the other term in (3.25) can be estimated similarly. By the definition of Q\ in (3.4), we 
have 


h = / ^ 

J 3 v v 

= I\ + if + If. 


2 

)uiy + Ji + ^= Y ]dy - J v\ 3 biRidy - e J vX 3 h ^ J £,f&{d£dy 

3 = 1 ' 


Since 




< C{5 + X 0 ){K\ + ||^f) + CSe 2 ( 1 + t)~ 1 E l + C{6 + X 0 )\\^i y \\ 2 + C6e b ( 1 +1)“5, 


2 , 


and 


J (| Ji| + |||) • |6i|dy < C{8 + A 0 )(/ii + ||$J 2 ) + C5e 2 ( 1 + t)" 1 ^ + C<fe 5 (l +*)"§, 


we obtain 


I\ < C{5 + A 0 )(iFr + ||<M| 2 ) + C(5 + A„)||tMi 2 
+ C6e?(l + t)~ 1 E 1 + C5e 2 ( 1 + 1)" 1 . 


On the other hand, from (2.37), we have 


(3.26) 


i?i = 0(l)5e 2 (l + t) x e 4c (®±H 1 + t ), as x —> ±oo. 


From (2.18), Q y satisfies 


\0 y \ = 0(8)e(l + t) 2 e 4 “( 9 ±)(i+*), as x —> ±oo. 


Thus, by (2.4) and the assumption on the profile, we have 


«( 0 ±) = \^( e ±) > \^±)- 


(3.27) 


Since a[9±) = 3 f ^g ^ , b(9±) = max{a(0±), } and c(9±) = max{a(0±), \b(9±)} , it follows 

from (3.27) that a(9±) > §c(0±), which leads to 


\lf\ < 


^ J \9 y \b 2 dy + CSe 2 (l + t)- 1 . 


(3.28) 
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We now estimate the integral if. Let M* be a global Maxwellian with the stat e 
satisfying \9 < 9* < 9 and \p — p*\ + \eu — u*| + \9 — 0*| < rjo so that Lemma 
that, 


3.3 


holds. Note 


/f =-e J v\ 3 bi j £j®\d£dy - e J vX 3 b 3 j £f®jd£dy =: if 1 + if 2 . (3.29) 

The estimation on if 1 is straightforward by using the intrinsic dissipation on b\ and (2.24[). 


If 1 \ = \-e I vX 3 h I tiLjti-PLfaGJ-eQiGMdtdyl 


— C J 1^1 I (j { £U yi Gy) | + \(£Uyy,6yy\ + \(£Uy,9y)\\(Vy,£Uy,9y)\Jdy 

<C5J \9y\bldy + C8e 2 {l + t)- 1 + C5\\{(j) y ^yXy)\\ 2 - (3.30) 

The estimation on if 2 is more complicated and it will be divided into five parts as follows. From 
(2.24), it holds that 

If 2 = -e J vX 3 b\ j gLj}(G r )d£dy + e J vX 3 h j&^Lj}(G y )d£dy 


-e / vX 3 b\- / eiL~ M [Pi{^Gy)]didy + e 2 / vX 3 b 3 / ^L^[Q(G,GMdy 


(3.31) 


+ 2e 2 / vX 3 b 3 / eiLlf[Q(G,GMdy =: 

J ^ i=i 

For the integral if 21 , one has 

j vX 3 h J gLj}(G T )d£dy - £ j vX 3 b\ f ^Lj}{G r )d^dy =: if 211 + if 212 . (3.32) 


If 21 = -e 


Note that the linearized operator satisfies that, for any h G TV 2 -, 

(L M l h) T = L~ M \h T ) - 2 Lj}{Q(Lj}h,M r )}, 
( I'm h)y = I'm ( h y ) — 2 L m {Q(L m h, M y )}. 


(3.33) 


Then it follows that 


lf 2U = -£ 


vX 3 bi J if{L^G) T didy - 2e J vX 3 bi J g Lj{Q{LjJG,M T )}d£dy 
= -(e J vX 3 b\ J ^L^fGd^dy) T + e J (uA 3 6i) t J £? L^Gd^dy (3.34) 


-2e / vX 3 bi / L M {Q(L M G,M T )}d£dy. 


The Holder inequality and Lemma 3.4 yield that 

| [ £lj}G<% I 2 < C [ eMlei) _1 M*de- f U -W\L-^G\ 2 dt < C 


Kieo,^ 


AT 


AT 


|G|X- 


Moreover, from Lemmas 3.2 3.4 one has 


gL-JiQiL-fG^Qjdtf <C [ ^\L^{Q(L^G,M T )}\ 2 d^ 


< C 


Kiel) 


-i 


Af< 

< C(u 2 + e 2 u 2 + 0 2 ) 


j M, '~ M 
\Q(L^G,M T )\ 2 dC<C 

r Kiei j- 1 


Kiei) lr -i/„2, t Miei).,,* 


AT 




AT 


|A/ r | 2 d£ (3.35) 


AT 


-|G|X- 
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Combining (3.34)-(3.35) gives that 


If 211 < - 


e [ vXsh [ tfL-jGdtdy) + cp\\ ($ T , tf T , W T )W 2 + CSe 2 (l + t)~ 1 E 1 

1 J Kiel),- - - (a36) 


+ C/3 £ 


M* 


-\Gfd%dy + C'AoK^tjV’tjC 


2 

ry || > 


where and in the sequel /3 is small positive constant to be chosen later and C /3 is a positive 
constant depending on (3. By the definition of G in (2.24), similar to the estimate in (3.30), one 
has 


lx I = I ej v\ 3 h J ttL- M \G T )dtdy\ 

< CS£ 2 (1 +t)~ 1 E 3 +Cd£ 3 (l+t)-l + G<5|| RRKr 


(3.37) 


Substituting (|3.36|) and (|3.37|) into (|3.32|) implies that 
I? 21 


‘i <- 


e J vX 3 h j ilL-^Gd^dy) r + C/3|| W T ) || 2 + C5£ 2 {1+ t)~ 1 E l 


+ C/3 £ 


KKP 

XL 


\G\ 2 d£dy + C(6 + A 0 )||(<K, W, Cr)|| 2 + C5£ 3 ( 1 + t)-2. 


(3.38) 


The estimation on J 32 * (i = 2,4,5) is straightforward by using the Cauchy inequality and 

KICI), 


Lemmas 3.2 3.4 First, it holds that 

Ilf 2 ! < CS£ 2 ( 1 + t)~ 1 E 1 + CXqK\ + C£ 2 


AL 


-I G y \ 2 d£dy 


+ Cfe 3 (l + f) 2 -f Cd£ 2 \\((j), i/j, Qy 


(3.39) 


Since 


£ 2 L-/{Q(G,G)R| 2 < C I ^®|L^{g(G,G)}| 2 de 


<c / ^Ml\Q(G,G)\ 2 d^<C f ^\L^G\ 2 dL 


M* 

<C\(£u x ,0 x )\ 2 


AL 




|G|R, 


and 


£ 2 L^{Q(G,G)R| <G( 


Kiei) lr _i 


M* 


L^{Q(G,G)}\ 2 dfr 


<C( 


Kiel) 


-l 


AL 


-|Q(g,g)| 2 #)* <<? 


Kiel) 

AG 


iG| 2 de, 


it follows that 

I/? 24 ! + I/? 25 ! < C{5 + A 0 )e 2 


Kiel) , ru-2/ 


-\G\ 2 didy + C5£ 2 (l + t)~ 1 E 1 . 


AL 


(3.40) 


The estimate on 1 323 is similar to the one for I 321 . First, notice that 

4 

PliHlGy) = {PifaG)}y + J2faG,Xj)Pl(Xjy)- 

3=0 
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From (3.33) and Lemmas 3.2|3.4 we have 


r323 _ 
~ 


e J (-A 3 6i) tf y gLjtfrfaGMdy-e J -Xsh J £fL^^^G^fPfXjyMdy 

j =o 

-2 e f -A 3 h [ {Q(L]^[Pi(£iG)\, M y )}d£dy 


<C« e 2 


Kiel), 


M* 


|G| 2 d£A/ + C<5e 2 (l + t) _1 £i + C(X 0 + P){I<i + ll^f) + CA 0 ||(^, $ y , ( y )\\ 2 , 

(3.41) 


where we have used the fact that 


K 6 e,*)|>< c/^«. 


Substituting (3.38), ( 3.39| ), ( 3.40[ ) and (3.41) into ( |3.31 ) gives that 

- (e f v\ 3 h j &L-£Gd£dy) + CSe 2 ( 1 + t)~ 1 E 1 + C(A 0 + P){K 4 + ||$J 2 ) 

Kiel) 


if < -( 


^^|G| 2 d£dy + Ce 2 


M* 


+ C/3|| (<h,4h!F) r || 2 + C>K 

+ C(5 + A 0 ) ^ IKKKKe)l | 2 + ^e 3 (l + t)' 
l«l=i 


M* 


|G y | 2 d£dy 


which implies by (|3.29|) and (|3.30|) that 
A 3 < - 


3 <-(ej vX 3 bi j gLjlGdtdy^ + C5e 2 (l + t)~ 1 E 1 + C(A 0 + P){K X + H^f) 

^\G y \ 2 dHdy (3.42) 


+ C'/3||(4>, 4', 1F) t || + C /3 e 


^M\G\ 2 dtdy + Ce 2 


+c(s+ A 0 ) ^ iid“(0,v>,e)ii 2 +c^ 2 (i+tr i . 

i“i=i 


And finally, (|3.26|), (|3.28|) and (|3.42|) yield the estimate on /, as follows. 

A < - | 


J vX 3 b\ J ffl/jfGd^dy^J + C<5e 2 (l + f) 1 Ai + C(Aq + + ||$y|| ) 


+ C/3||(4 > ,4/,1F ) t || 2 + C / j e 2 


Kiei),^2 JtJ „., ^2 


M* 


Kiel), 


|G| 2 ^A/ + Ce 2 / / ^|G y | 2 d^ ( 3 . 43 ) 


+ C(5 + A 0 ) £ ||5“(0,^e)H 2 + A<5e 2 (l + t )- 1 + ^ J \6 y \b 2 dy , 


l«l=i 


which completes the estimate on the term / vQi^idx. 


Estimate on 


q 2 b 2 dy: 


Notice that the profile has no intrinsic dissipation on b 2 . Fortunately, it holds that q 2 = 
and Q 4 has the decay rate as e 3 (l + t)~ §. Thus the estimation on f q 2 b 2 dy can be 


5 p+ 

directly obtained even though there is no intrinsic dissipation on b 2 . For example, 
| [ £U\Rib 2 dy\ < G<5e 2 (l + tf l E\ + Cfe 3 (l + t)~f 

e 2 u 1 b 2 e i Q 1 dHdy\<C6£ 2 (l + t)- 1 E 1 +C(5 + X 0 ) £ ||d“(</>, C )|| 2 

|a| = l 


+ C5e 1 


Klfl )G* 

EL 


d£dy + C 


KKP. 

EL 


(G 2 + G 2 )d^dy + G5e 3 ( 1 + t) 2 . 































25 


And the term e / / £i|£| 2 0 f 62 d£dy can be estimated similarly as for if 2 where the intrinsic 


dissipation on b\, 63 is not needed. Notice also that all the other terms in q<± are of higher order. 
Therefore, one has 


I 2 — j I'lb'zdldy < (e J J A(£,b 2 )L f jGd£ : dy) T + C5e 2 (1 + t) l E\ + C*(Ao + @)(Ki + ||4>, / || 2 ) 
+ C/3\\($,'S>,W) T \\ 2 + Cpe 2 f J"M\G\^dy + Ce>J J ^\G y \ 2 d^dy 
+ C(5 + A 0 ) Y \\d a {(t>^,Q\\ 2 + C5e 2 (l + t)~\ (3.44) 

|a| = l 


where A(£, 62 ) is a linear function of 62 and a polynomial function of £. Using (3.43), (3.44) and 
(3.21), we get 


Eir+(J J eii(£, S)L^Gdedj/) T + + J \0 y \{b\ + 6|)dy 

< Cfe 2 (l + t)~ l Ex + C(A 0 + /3)(A'i + ||$J 2 ) + C/?||($, , W) T \\ 2 

+ Cpe 2 J j^^didy + Ce 2 j J u{ ^\G y \ 2 d^dy 

+ C(5 + A 0 ) Y \\d a {<t>,il>X)\\ 2 + C5e 2 (l + t)-\ 

\a\=l 


(3.45) 


where we have used the smallness of 5 and Eo- Here A\ is a linear function of B = ( 61 , 62 , 63 )* 
and a polynomial function of £. 

Note that K\ does not contain the norm || < h y || 2 . To complete the lower order inequality, we 


have to estimate 4?^. From (3.3) 2 , we have 


- T 1t + = A Wy - ^%£ 2 N ly - Q V 

3v y v 3v u 15 p+v y 


(3.46) 


Multiplying (3.46) by yields 


2 d(d) ^2^ / 2 M 0 )\ 2 , P+k2 _( 2 t it 8 d(d) 2 


( ^ *»> T ~ ~ + = ( 


3v y 15p+u 


Nxy - QlJ 




Since 


“ ($r^l)i, + ^ly ~ -= -e^l^i y, 

op+ 


we can obtain 


(/^$ 2 -^!^) r + J ^ 2 y dy <C\\^ ly ,W y )\\ 2 + Cde 3 (l + t)- 3/2 + J Q\dy. (3.47) 


The formula (3.4) for Q\ and the Cauchy inequality directly yield 


g?^<C(6 + Ao)(AA + ||«h,|| 2 ) + C'A 0 Y 

\a\=l 

+ Cde 3 ( 1 + t)~ 3/2 + C [\ [ gQxdtfdy. 


(3.48) 
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And using Lemmas 3.2 3.4 implies 


gQidtfdtdyKCe 2 J J ^M(\G y \ 2 + \G T \ 2 )d^dy + C (6 + X 0 )e A J f <%dy 


+ C(S + A 0 ) Y \\d a (<j), if), C)|| 2 + CSe 3 (l + t)~V 2 . 

\a\=l 


(3.49) 


Substituting ( |3.48 ) and (3.49) into (3.47) yields 

<c 2 k 1 +c 2 £ 2 r 


„ (| G y \ 2 + \G T \ 2 )dtdy + C 2 (5 + \ 0 )e 

V «/ dVL* „ 

+ c 2 ( 5 + a 0 ) y ii9“(0,^c)ii 2 +^ 3 (i+r 3/2 - 

|o|=l 


'(lfl)|Gf 

M* 


d£dy (3.50) 


Multiplying (3.5) by e 2 ](f-, one can obtain 


, 2 . 2 vG ~ f 1 n r/- / I? — ™| 2 1 > 1 

(£ 2 M? t ~ £ M* LmG ~\R 9 Pl[il{ 29 + 


,)M] 


+ £U\Gy — -Pi(£iGy) + £vQ(G , G) — uG t |’ • 


6' 

M* 


(3.51) 


Integrating (3.51) with respect to £ and y and using the Cauchy inequality and Lemmas 3.2 3.4 
one has 


1 


2 M, 


-|G| 2 d£dy) 


3cr 2 

+ T e 

4 


Kiei),/v, a 


M* 


|G| 2 d£dy 


<C3<5e 3 (l + t)- 3/2 + G 3 £ ||<9"(^>, if), C)|| 2 + G 3 £' 

\a\=l 


K|g|) 

Af* 


|G y | 2 d£dy. 


(3.52) 


On the other hand, since (4>,VL) r can be represented by (4 > ,4',bL)y and (4>, dL W) yy from 
the equation (3.3), we can get an estimate for (4>, VL) T as follows. 


\\^,^,W) T \\ 2 <C 4 (K 1 + \\^ y \\ 2 ) + C 4 Y \\d a {<t>AA )\\ 2 + C i 5£\l + t)~ 

l«l=i 


+ C 4 ^ 


Kl^l) 

Af* 


(|G y | 2 + |G r | 2 )d£dy + 04(6 + Aq)K 


Kiei)|gp 

Af* 


■d£dy. (3.53) 


Now we can complete the lower order estimate. Since Ai is a linear function of the vector B 
and a polynomial of £, we get 


A^B^Gdidyl < -E\ + C £ 2 


\Gf 

A L 


d£dy. 


We choose large constants C 3 > 1, C 2 > 1, C 3 > 1 and small constant f3 such that 


CiEi+Cxe Ij AiL^Gdidy + C 2 J (^U 2 - %^i)dy + C 3 e 

>\c>E I + C 2 l + ^ f [ 


\Gf 

2M, 


-d£dy 
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C 2 C 2 -C l C i C 4 )K l + 


(I- 


s T K ' + <?2 /S*>' 


C 2 ^-(5 + /3 + A 0 )CiCi(1 + C 4 ) 

4t; 




and 


|C 3 - CiCiC 4 (<5 + p + A 0 ) - CpC^ - + A 0 ) > jC 3 . 


Hence, by multiplying (3.45) by C\, (3.50) by C 2 , (3.52) by C 3 , (3.53) by C\(5 + Eq + E\)C\ and 
adding all these inequalities together, we have 


E 2t + Aa + J \6 y \(b 2 b 2 )dy < C§5 e 2 (11) 1 (-E'a4-1) 
+ C§£ 


_ 2 / / Kiel) 


M* 


where 


(\G y \ 2 + \G T \ 2 )dHdy + C 5 \\d a {^X)\\\ (3.54) 

l“l=i 


|G | 2 


E 2 = C\E\ + Ci J J EA^Gdidy + C 2 J - QyV^dy + C 3 e 2 J J ^-d£,dy , 


K 2 = ( ±K 1 + C 2 j ^L$ldy+\\(<S>^,W) T \\ 2 + ^C 3 £ 2 


KieDiKp 

M* 


d£dy. 


3.3 Derivative estimate 

In this subsection, we derive the higher order estimate for (<3?, 4', W). First, by the definition of 


©J in (2.24), it holds that 

-e J = £ 2 IV 4 + e 3 F 4 , 

-^4 = f41@x@x H - f^Vx^x H - f^S^x ~F f44@xxi 

J-RjJ = ^(1)[(0#| \0x\ "F | @x | \eu x \ | SU X |) | Ux | “1“ |^a;||^x| “1“ |^®a;|]* 

From (2.9) and ( |2.36| ), one has 


(3.55) 


On" 

0r “01 y = -^1 yi 

4 e ai{ 6 ) fi( 0 ) _ 

01T H“ (p P)y q ( - 'U'lyjy + Q 5 > 

o V 1 ) 

Ipir — £( 'Uiy Z ^iy)y 4“ ^4+z? 2 — 2, 3, 

V V 

Cr 4“ EpU]y Epit\y — (— Oy ~zz $y )y 4“ ^8? 


(3.56) 


where 


Q 5 — £ 2 (-/V 4 J ! — N 4 y) + £ 3 F 4 y — £ J gelyd^ 4“ (£~N 4 y — R\y) , 
Q 4+ i = £ 2 (AK - iVjy) 4- £ 3 F iy -£ £i£i@lyd4 - Riy, * = 2,3, 


q 8 = h 2 Jv„+ e\n 1v - N ly )+ s 3 f 1s + i -d£l £Vlv + ^2 fEUEu 


p( 0 ) 


3 v 


*y 


0 f Kiel 2 ©nKe+E 

J i=i 


i =2 


£Ui eiKQly^e -Hy -{- ^(l^^l )r 4“ EPyU l, 
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and Ni,Fi (i = 1,2, 3,4) are defined in (2.26), ( 2.31| ) and (3.55) respectively and N t (z = 1,2,3,4) 
is the corresponding function of IVj (i = 1,2, 3,4) by substituting the variable (v,u,6) by the 
profile ( v , u, 0). 

We will use the convex entropy for the fluid system to obtain the first-order derivative 


estimates of ($ y , 4> y , W y ). Multiplying (3.56) 2 by ipi and (3.56) s by tfti, one has 


,1 V-^ , 9 \ , . 4e / /i(0) _ \ ,/r(0) /r(0) - \ , ri , / 

( 7 ^ / Vi Jr \P P)Vly 4 “ ^ ~ ^1 y ~z ^ly J Vl y 4 " ^ Z H iy)Viy — / Qa+iyi 4 “ (' 

2=1 1=1 


Since p-p = - |) 4- we obtain 


, 1 v—^ 2 \ 2-1 l w 2 4/i(0) )2 s^p(6) ,2 4e,/i(0) /x(0). _ , 

Z^ ly + 3 v ^ + „ ^h/ + 3 ( „ ,u J^iyV'i 


'2 fr' 3 v 

,p(6) p(6) 


i=l 

3 


3n 

3 


i=2 


3 V 




_^ ^ i 

)Uiylpiy = ^ ' p’iQi+A 4 “0( — z J^Vij/ 4“ (' ' ' ) 


i=2 


«=1 


5 u tb 


(3.57) 


Let 


<f>(s) = s — 1 — Ins, 


then it holds that 

{^<f>(z)}r = \oM”) + \e{-~ + z)0r + \K~ 2 + r)u* + ^(-- 4- zK 
3 v 3 v 3 v v 3 v z v 3 v v 

= -0(-- + \)(t>T ~P^C^)vt + upr^(z), (3.58) 

3 v v v v 


where 


'L(s) = s 1 — 1 + Ins. 

It is easy to check that 4>(1) = 6'(1) = 41(1) = \L'(1) = 0 and 4>(s), 'L(s) are strictly convex 


around s = 1. Substituting (3.58) into (3.57) yields that 

VV + lm- 


1 , 0 2 ,v s \ 2 „ , 4 p(6) 2 2 4e u(0) /x(0) 


2=1 

3 


-^ + i^ + £^4 + T( 


3 u 


1=2 


3 v v 




+£Y'(— - ^^-)u iy tpi y = V i>iQi+A 4- ^-0(- - z)-/Viy -/^(zjw 4- fPrd>(4). (3.59) 

' 7 J 7 ? * * ^ h 7 J 7 J * 7 J 7 J 


2 = 2 


v v ’ " " z —' ' 5 v v 

On the other hand, multiplying (3.56 ) 4 by (, it holds that 

c 


C _(<0) o K (0) s ^ C +Qg c_ 


^Cr + e(puiy - pui„)| = (^9y - ^y) 


v y 7y0 "°e 


(3.60) 


One can compute that 
Ur = 


+ 0 r ^(!) = (04»(=))^ + O(l)fe 2 (l + f)- 1 |C| 2 , (3.61) 


s(pu ly - pu ly )-jj = ^ipiy + e(p - p)u ly ^ = ^i/j ly + 0(l)5£ 2 (l + t) ^(V’jOI 2 ! (3.62) 
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and 


K (°) p _ h ^ln\ C , N 6>At((9) 2 «(g)gyCCy MyCy ~ \°y\ 2 ( ( K (°) K (°) N 

u y V y JyO y v6 2 y v6 2 $ 2 ' V V 


< 


3fe W c j + Cfc 2 ( l + tJ-'KAOl 2 . 


^ 4 ufl 2 

Substituting (3.61), ( 3.62| ) and (3.63) into (3.60) yields that 

K>) 


, 2 C , , ^ 

+ 3^ ly+ 4 u0 2 ^ 


^ \ y < (• • • )y + C5e 2 (l + t) 1 \(4>, C)| 2 + |Q8^|- 


Combining (3.64) and (3.59) and using Cauchy inequality, one has 
E 3t + \l < 3 < C5s\ 1 + t)- 1 ^ + C5e\ 1 + f)“ 2 + Ce 2 (l + f)" 1 


l4|(6 2 + 5 2 )dy 


+ 


K^-Wll/l - — 
vv 


<t>+^)dy +^2 l J *+2|> 


2 — 1 


where 




k 3=i(i *bi, +s* +1 : m 2 ) <&. 


3 ^W,/,2 , 3 g*(fl) /-2 


and 


h= £ 2 (N 4y - N±y)ipidy + / E 2 F iy ilj l dy - £ 




^2+i — 


J e 2 (N iy 



- N t y)'iJj t dy + / e 6 Fiy^dy - e / / €i£i@i y ipid£dy, i = 2,3, 


,C 


C 


h= E ( N iy - Niy)-dy + / e FWdy - -£ 




+ E 


£ 2 u?; 


Z=1 




(3.63) 


(3.64) 


(3.65) 


(3.66) 


In the estimate of (3.65), we have used the estimate like 

| J(E 2 N 4y - R ly )il)idy\ = | J(E 2 N iy - Ri y )^iydy\ < C J \(E 2 N 4yy - ^i^KI&iI + N)<fr/ 

< C5e^( 1 + t) 2 + C£ 2 (l + t) f \d y \(b\ + b 2 )dy. 


Now, we calculate the terms on the right hand side of (3.65). Firstly, a direct calculation yields 

— ~^z4> + ^ — & y ) + 0(1) |(0, 4>, C)| 2 + le^ii/^il + <fe(l + t) 1 


and thus by combining (3.23) and (3.24), it holds that 

2 .2 aV / 


5- -'M-^+l)* 1 


-£ IVi„ - — 


< Cfe (1 +1)- 1 ^ + C5e\ 1 + f)" 2 + CI e 2 [\(N ly --) y \ + |JVi w e«i y |] • |(6i, b 3 )\dy 


< C5e 2 ( 1 +1)-^ 3 + Cfe 4 (l + t)- 2 + Ce 2 (l + t)" 1 / \9 y \(b{ + b 2 )dy 


(3.67) 
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By using the definition of N. 4 and F 4 in (3.55), one can obtain 


I £ 2 (N 4: 


1 


y - N iy )^idy + / e i F iy ^idy\ < —if 3 + / e 4 (iV 4 - N 4 ) z dy + e b F£dy 


1 


< ^K 3 + CdUyf + C5s 2 ( 1 +t)" 1 ^ + C 6 e\l + t)" 2 . 


(3.68) 


And by using Lemma 3.2 3.4 one has that 


< —if 3 + Ce 2 
“32 J 


ei®ldt\ 2 dy 


< —if 3 + Ce 2 
“32 6 


Kiel) 


M, 

4/i , .1-2 , r<A 


\{G y ,G T )\ 2 didy + C5 Y \\d a (cj),iP,0 \\ 2 + C 6 e 2 (l+t)- 1 E 3 
l«|=l 


+ C5e 4 ( 1 +f)" 2 + Ce 4 5(1 +f)" 1 + 


\Gf 

M* 


den 


L°° 


Kiei)|g| s 

M* 




(3.69) 


Combining (3.68) and (3.69) yields that 


K < 32^ 3 ^ e2 


/ / ^pKG^GOl'dedy + Cd ^ H^^Of + ^d + t)- 1 ^ 

J J * M=i 


+ C5e 4 (l + t)~ 2 + Ce 4 <J(1 + f)" 1 + 


|K | 2 


M, 


Klh 


KieDlKp 

M* 


d£dy. (3.70) 


Similarly, I 4 .I $. Iq can be controlled by the right hand side of (3.70). Substituting (3.67) and 


(3.70) into (3.65) gives that 


e 3t + 2^3 < 


C 6 s 2 [ j ^\(G y ,G T )\ 2 dtdy + C 6 6 Y IKKKKOII 2 


l«l=l 

+ C6<fe 2 (1 + f) 4 E 3 + Cq 8 £^{\ + 1 ) + Cq £ 2 (1 + 1 ) 4 f \9y\ (6| + b 2 )dy 


+ Cq£ 4 


5(1 + 1 ) + 


\&t 

M* 


deih 


Kiei)ici 5 

Ad* 


d£dy. 


(3.71) 


Note that the norm IK^H is not included in K 3 (see ( 3.66[ )). To complete the first-order 
derivative estimate, we follow the same way as to estimate in the previous section. By using 


the (3.56 ) 1 , we can rewrite the equation (3.56 ) 2 as 


4 p(0) 4 e 2 y(0) ~ 4 m{ 6 ) 

d>yr — ipir ~ (P ~ P)y — ^ z—Niyy — — ( — ) y il>i y 


3 v 


3 v 


f i^ Uly ~ ^V^ lv) \ y +£ i + R lV ( 3 - 72 ) 


Multiplying the equation (3.72) by <p v , one has 


2 n( 0 ) 2 2 n( 6 ) 2 , , , _x , 

q ( < t > y) T ~ O ( = )r4 > y ~ ~ (.P ~ P)y4 > y 


3 U 


3 v 

r 4 s 2 y( 6 ) ~ A ai(9) 4e w p(9) p(9) / 2 - i 

j Z—Niyy — -(—— )y^ly ^“[( — ~ £—) u ly\y + e / £l@lyd£ + R\ y j(f)y. 


(3.73) 


/ -N P, 2 p p 2 1 1 

-{P-P)y = ~<Py ~ w=(y + K ~ 3)% - ,(- - 1 ^ s , 
u 3u v v 3 v v 


Since 
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and 


3e 2 


(frytplT = (4>y1pl )r ~ + 1p( y - — N ly ^ ly , 

o 

integrating (3.73|) with respect to y and using the Cauchy inequality yield 


2 m 


P j.2 


+ I -Zz&dy < C 7 K 3 + C 7 8e\ 1 + + C 7 5e\ 1 + t)"5 


27 


37 ^ 

+ C 7 (S + A 0 )e X] \\d a ^X)\\ 2 + C 7 5\\d T ^,Q\\ 2 + Ce 2 j \ j&e ly dt\ 2 dy, (3.74) 


M=2 

where we have used the fact that 




V V 


y w + C < 5 e 2 (l + t) 1 £3 + C A 3 . 


It follows from (2.24) and Lemmas 3.2 3.3 that 



ei®i y di\ 2 dy<C5e b (l + t)-i+C5e\l + t)- 1 ^ ||d a (^,C)l| S 

|o|=l 

+ Cfe 2 £ \\d a (cj>, ip, C)l| 2 + CJ 3 , 

\ a \=2 


where 

Js = 


■ |«|=2 


K|g|) 

M* 


\d a G\ 2 dtdy + e\5 + \ 0 ) £ [ f ^\d a G\ 2 d^dy 

\a\=l J J * 


+ ^ + () -. + // ^ + / / %£«,) / / 


To estimate (0, ip,() T , we use (3.56) to obtain 


R(0, ip, C)|| 2 < C 8 (£ 3 + H^ll 2 ) + C 8 <5e 2 (l + f)- 1 ^ + COe 5 (l +1)"§ 
+ c 8 X] ||9 q (0^,0II 2 + c 8 j 3 . 


I«l=2 


Thus we choose large constants C 4 and C 5 so that 

C 4 E 3 + C 5 J ( 2l ^ ] 0 2 tVO) dy > £.3 + C 5 J ^rfdy, 


and 


Let 


- C5C7 - C 8 > Jc 4 , 0 5 [ £fidy - C 8 ||^|| 2 > ^ [ l<gdy. 


27 


4 / v 


£4 — 646 ^£3 + C*5£ 


.-2 




(3.75) 


(3.76) 


(3.77) 


£4 = *C 4£“ 2 £'3 + ^e- 2 [ ^cftdy + £ —2 


> Vv> Cr) || 2 - 
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Then from ( pT7l| ), d3l74| >, ( ]375j ) and ([3777]), we have the following estimate on the (<p, ip, 0 


E 4T + K A <C 9 6£ 2 (l + t)~ 1 E 4 + C 9 S£ 2 {l + t)- 2 + C 9 e- 2 Y ||<9 Q (0 iP, C)lP 

M=2 

+ C'g(l+t) 1 f \9y\ (b 2 + b^)dy + Cg£ 2 J 3 , 


(3.78) 


where J 3 is defined in (3.76). 
Define 

E-j = E 4 + £ 


m^dy, K, = K, + °-e 


-(If l),A 2 


M* 


\G\ 2 d£dy. 


(3.79) 


Then from (3.78) and ( 3.52| ), one has 


E bT + Ks<C l9 5£ 2 (l + t)- l E b + C l0 5£ 2 (l + t)- 3 2+C W £~ 2 Y \\d a ((P,ip,0\\ 2 

\ a \=2 

+ C w (5 + X 0 )£ 2 Y I j ^ \d a G\ 2 d£dy + C w Y J J ^\d a G\ 2 d£dy 


|o| = l • 


+ Cio(l + 0 1 J \0y\(Pi + b 2 )dy. 


I«l= 2 ‘ 


Next we derive the higher order derivative estimate. Applying d y to (3.56) yields that 


where 


3 2 - 

(j)yr Vhl/J/ - El\yy, 

/ . P+/- P+A _ Ae fp( e ) MW- X , n 

Wlyr + C yy - Pyy ~ o ( ,, Ul V T . u ly)yy + 

u u o u u 

'Ipiyr — £•( “ 'U'iy Z ^ iy)yy ^ — 2, 3, 

C yr + P+'tplyy = (“—~@y - -pr^dy)y + ^^Eiy + Q\2, 

V V 0 


Q g ~ y +( ^yy + (~ ~ ~)^yy + G(i)(\v yy \ ■ \(<p, 01 + 1 4*Cyy I ) 


(3.80) 


3 ( v 2 v v e y -- 


1 _ , \ 4 / 0 2 1 2 

- ,r " ") ' :>,( 


o\ v 3 V y Q3 V y) e I 0 G>\yydE, Rlyy, 


Qi+ 8 — £ P.l^i^dlyydP, Ri yy ,i — 2,3, 

Ql 2 = ~£Utyy(p ~ p) - £{PyUy ~ pylly ) + (j)+ ~ p)lplyy + Q\ 3 y 

+ Y £2 ( Ui J y d 0 


1 

~ 2 £ 


2—1 


4 MO 22 

3 u ly 


Qi 3 = x — e~u 


+ e ' 2 E 


i=2 


^p~ u ly - H iy - R 4y + \{\ £u \\ + epj/Ui- 
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Multiplying (3.80) 1 by p+4> y , (3.80) 2 by vip\ y , (3.80) 3 by ipi y , (3.80) 4 by ( y , we have 


+ 2^ y + ^ ^ + 9.^ dy 


i=2 


j t 4 / L 3u 


i=2 


dy 


<C5e(l + t) 2 \\(<f>y ,1py, Cy)\\ 2 + C5e 3 (l + t) 2 IK^) VOII 2 + C\\(4>, 


yi Vfyi Cy 


|£| 3 eiyd£ dy 


< 


+ C|| ip y , 0)11 3 + C(5e' 6 (l +t) 3 + Ce 2 j 

C6e(l + t)~^ ^2 jj|d a (000ll 2 + Ctfe 3 (l + t)“§ ||(0, -0)|| 2 + C\\ (0,-0)|| ||(<0> '•PyXy) II 3 

|a|=l 

+ C'||(^,^,Cy)ll^+Cfe 5 (l + t)-i+Cfe 2 £ ||d“(</>,V’,C)l| 2 + C'J3, (3.81) 

H=2 


where we have used (3.75) in the last inequality. 
Let 


E 6 = 


P+ ,2 , v / 2 , ,2 , ■ L > 

<0 + ^ + 2 ^ + ^ 


i 


i=2 


2 V i 


dy, Kg = 


i=2 


3v 


dy, 


then (3.81) implies 


E 6t + \k % < C u 5e(l + t)~* J2 O Of + Cn<5e 3 (l + 1)~ § iP, Of + C n Se 5 ( 1 + t)“ 

M =1 

+ C 1 i||(</-,V>)||||(^,^,C0l | 3 + Cii||(^,^,C0ll^+Cnfe 2 2 IO“(OOOH 2 + CiiJ 3 . (3.82) 

M=2 


To get the estimate on <j> yy , we use the momentum equation (2.8 ) 2 - Applying d y on (2.8) 


2 ) 


it holds that 


Iplyr T (p P)yy T £U]_yr T Pyy — £ 0 Gyyd^. 


(3.83) 


Note that 


/ p , 2 1 . _ 4> _ 2v v . 2p v 

(P - P)yy = --<t>yy + ^ yy ~v^ P ~ P ’ Vyy ~ “ ~V^ ~ p ’ y ~ V^’ 


then multiplying (3.83) by — cp yy and integrating the reduced equation with respect to y give 
that 


P I 2 7 


( j *l ! ,y4 > yydy)T 

< 6*12^6 + Ci2<5e(l + t) 2 11(00 C)j/1| 2 + C'i2(5e 3 (l + i) 2 IKOOOf T Ci2<fe 5 (l + i) 5 

+ C'i2||(^,^,0)ll^+C'i2e 2 f t U -W\G yy \ 2 dtdy. (3.84) 


M* 


To estimate (0 ip, C,) yT and (<; t>, ip, () TT , we also use the original fluid-type equation (2.8). Here we 
only consider the term ^ ip“l yT dy because the other terms can be estimated similarly. It follows 
from ( 2 . 8) 2 that 


Iplyr — C P p)yy ^lllyr Pyy & / £>\Gyyd£,- 


(3.85) 


to|cn 
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By (3.85) and using the Cauchy inequality, it holds that 


UiyrW 2 < C 13 (K 6 + UyyW 2 ) + C 13 Se( 1 + t )-2 II (</>, y ,, Oy II 2 + c 13 5e 3 ( 1 + t)~* || (</>, </>, C)||' : 


+ C'i3^£ 5 (l + t) 2 + C 13 II ((f>y, l/jy, Cy)II 3 + Cl 3 £ 


Kiel) 

M, 10yy 


\G yy \ 2 dtdy. (3.86) 


Let Cq and C 7 be suitably large constants, then it follows from (3.82), (3.84) and (3.86) that 
C 7 (c 6 E 6 - j ihytyydy)^ Y \\d a (c/>, 0f 


< 


Ci4ds(l-\-t) 2 'y ' ||r)“( 0 , ip, C) ||“ + Ci4<5e 3 (l + t) 2 1| (</>, 4 > i C) l | 2 4" Gu<5£ 5 (1 + t) 
M =1 

+ C *14 || ( 0 , -0)11 II (0J/S V’j/, Cy)l| 3 + Cull (</>y, l/jy, Cy) II 3 + C 14 J3, 


(3.87) 


where J 3 is defined in ( 3.76| ). 

To close the a priori argument, we need to estimate the non-fluid component d a G, |a| = 1 , 2 . 
Applying d y on (3.5), we have 

vGy T — vL M G y = —VyG T + v v LmG + 2 Q(M y , G) — |-^-Pi[£i(^ ———Cy + £ ■ ~'4 , y )^]| 


+ j^ewiGy — P\{^iGy) + £vQ{G ,G) — uG r . 


(3.88) 


Multiplying (3.88) by then integrating the reduced equation with respect to £ and y and 
using the Cauchy inequality and Lemmas |3.2||3.4[ we have 

v \&v\ 2 f f Kiel),/, 12 


2 AC 


| a |=2 


-d£dyj 

Kiel) 


+ T 


AC 


| Gy | 


M |d“G| 2 d£dy + G 3 5e 3 (l + 1)~ 5/2 + G 3 e -2 £ ||3“(<£,Oil" 

* l«l=2 

+ G 3 i 5(1 + t ) _1 £ ||0‘W,C)|| 2 + C 3 e - 2 £ ||O“(0,^e)H 6 


|a|=l 

+ g 3 Y lia>,M)lllo 
H =1 


AG 


|a| = l 


Similarly, we can obtain the estimate for G r . Hence, one obtains that 
u|3 Q G | 2 _ , \ 3d 


E 

|or[=l 


2 AT 

Kiel) 


«)cf E 

M=i 


Kiel),*,™ 


AC 


|d"G|XC/ 


<C Y I \ ^\d a G\ 2 dtdy + C 3 6e 3 (l + t)- 5 / 2 + C 3 s- 2 Y IKKKKC)I|- 

\ a \= 2 J J * H =2 

+ g 3 <5(i + 1) _1 Y IKKKKC)II 2 + c 3 £- 2 Y ikkkkoii 6 

l“l=l M=1 

Kiei), 


+ g 3 Y 

|o|=l 


AC 


-\Gfd£dy. 


(3.89) 
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r t* is ( 1^11 I* 

Finally, we need the highest order estimate to control E MP \d a G\ 2 d£dy and / V’i y <S>y y dy 

|a|=2 J ■> * J 

in (|ih87|). To estimate f ipiy^ydy, it is sufficient to study the a priori estimate for ^ 

J M=2 


v\d a f\ 

2 M* 


due to (3.12) and (3.13). Applying d a , |or| = 2 to (3.8), one obtains that 


vd a f T - EvL M d a G - Euid a f y + £id a f y =-d a vf T + £d a uif y - E [<9 Q ^vd^f T - sd a ^u\d^f y 

l/3|=i 

+ Ed a [vLmG - vL^Go] + E 2 d a [vQ{G , G) - vQ{G 0 , G 0 )] + d a [ - 4>f T + i/jf y - evR^ . (3.90) 

d a f . 


Multiplying (3.90) by integrating the reduced equation with respect to £ and y and using 

±V1* 

the Cauchy inequality and Lemmas 3.2||3.5[ similar to the argument used in [24] . one gets that 
v\d a f\ 2 ,_ \ 3 a 


E 

l«l=2 


2 M* 


- d ^y) r + T E £2 

l«l=2 


"(lei),*,™ 


M* 


| d a G\ z d£dy 


10 


<C 3 5E 5 (l + t)- 5 / 2 + C 3 (6 + Vo + X^J2 ||d“(<PtPC)ll 2 + c E !I<9 Q (<^,C)II^ 

M= 2 M =1 

+ C < 3 ( 5 e(l + t) 2 ^2 ll^ Q ( 0 ) Vb C)l | 2 + C 3 <fe 3 (l + t) 2 11 ( 0 , Ip, C )|| 2 

|a|=l 

2 




<*V(i + t)- 2 + E ||a"((/>,^,c)ir 

|o| = l 

KICD./V ,2 


Kiel) ,^,2 


M* 




+ G 3 (S + 770 + Ag )e 2 


M* 1 y 

Choose large constants Cs > 1 and C 9 > 1 such that 

CsC 7 /^ 6£/g _ ^ ipi y (p yy ayj -+- - 

^ M =1 




(3.91) 


Ej = 


(c 6 E 6 - J M „dy ) + ) E / / + 5 E 


M=2 


v| 0 p /| 2 

2 M* 




E 6 


(llWM^Oyf + E ll< 9 "(<PV’,C)|| 2 + X] / / 

H =2 H= 2 J J 


i«i=i 


MM2 

^“Gl 2 


2M* 


d£dy-C<5(l + t)-2. 


Let 


AM = 


4e 3 


E ii8“w.*oii a + ^ e 


W -2 


l<|a|<2 


KKP, 

M* 


|9“G| 2 d£dy. 


Then from (3.87), (3.89) and (3.91), one obtains that 


£ 7T + iL 7 <Cfe 2 (l + f)- 5 / 2 + c4 E 11^.^011^+C'<5(l + i)-i||(0,^C)ll 2 


|a|=l 


+ C 


+ C 


5(l + t) a + E|(^>,-0)11 - ||(0j/, "02/, C?/)lll E ill^C^^OII 


(3.92) 


M=i 


<fe(l + f) x + E ( _ ll^ a (^)"0)C)l| 2 + £ 


l“l=i 


|d“G | 2 


d£<A/) 


Kiel) 

M* 


|G| 2 d£dy. 
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From (3.78), (3.79) and (3.92) and using the smallness of 5, Aq and e, we have 


{E a + E 7 ) r + \{Ka + K 7 ) < C5e 2 { 1 + t)- l Ei + C5e 2 { 1 + 1)~ 2 + C( 1 + t )" 1 f \9 y \ (bj + b 2 3 )dy 


+ C 


and 


8e(l + t) 1 + ^2 (^\\d a {^,ilj,0\\ 2 + e J I ^ d£dy) 

|a| = l 


(E 5 + E 7 ) T + ~(K 5 + I< 7 ) 


» m )^2 


M* 


\G\*d£dy, (3.93) 


< C5e (1 + t) E§ + C5e (1 + t) 2 + C{ 1 + t) / |0y|(lq + b 3 )dy. (3.94) 


4 The Proof of Main Result 


For a suitable large constant Cg, by combining (3.54) and (3.94) and using the smallness of 5, Ao 
and e, we have 


where 


Note that 


Eg T + Kg < C5e (1 + e~t ) Eg + CSe (1 + £ t ) , 


Eg = C 9 E 2 + E 5 + Eh, Kg = -(IC 2 + I< 5 + K 7 ) + / 1^1(6? + b\)dy. 


Eg >\\(^,W)\\ 2 +{^1(^,0^ 


+ £ 


v\Gf 

M* 


d£dy} 


+ < % 


(ll(0,-0,Oyl! 2 + Y ll<9"(^V5C)ll 2 + Y j J^p~ d ^ dy ) 


+ 7^ 

|o| = l 


l«M |a|=2 

\d a G\ 2 


2 M* 


d£dy — C5(l + e t) 2 >, 


Kg> Y \\d y (<S>,*,W)\\ 2 + c 2 \^ E \\d a (<t>AA)\? + e I f Y^±\G\ 2 dt;dy 


=i 


M=i 


2M* 


+ (5Eii a ”(**c)ii 2 + i E 1 ■ 1 ^ 

[q| =2 |ck|=1,2 


M \d a Gfd£dy > + - / \9 y \(b{+b 3 )dy, 


and 


e 2 E 7 < C5e 2 (1 + e 2 t)- § + C{K 4 + K 7 ), and e 2 (E 5 + E 7 ) < C5e 2 { 1 + £ 2 r)"f + CKg. (4.1) 
Then the Gronwall inequality yields that 

Eg < CV~8{ 1 + e 2 t) C oV *, [ Kgds < CV8( 1 + e 2 t) C oV *. (4.2) 

Jo 

Hence, it holds that 

||($, T ,! F )|| 2 < CV8(1 + e 2 t) C oV *. ( 4 . 3 ) 
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Multiplying (3.94) by (1 + e 2 r) gives 


[(1 + £ 2 t)(E 5 + E 7 )] t + -(1 + £ 2 t)(A 5 + K 7 ) 

< Cs 2 (E 5 + E 7 ) + C5e 2 { 1 + e 2 r)-^ + C 9 [ \Q y \(b\ + b 2 3 )dy. 


Integrating (4.4) with respect to r and using (4.2) and (4.1), one has that 


which yields 
In particular, one has 


(1 + £~t)(E 3 + E 7 ) + f — (1 + e 2 s )(/\5 + K 7 )ds 
Jo z 

< Ce 2 [ (E 5 + E 7 )ds + CVS( 1 + e 2 r)5 

Jo 

< Cy/S(l + e 2 r )2 + C j K$ds < CVS(1 + e 2 r) 2 , 

Jo 

(E 5 + E 7 ) < C\/5(l + e 2 r)-i 


0 ! 

AA 


■d£dy < CVS(1 + e 2 r) 2. 


On the other hand, multiplying (3.94) by (1 + e 2 r) 2, it holds 


(1 + £ 2 s)^(K 5 + K 7 )ds < CV6( 1 + £ 2 t) C ov/ 1 


Multiplying (3.93) by (1 + e 2 r) and using (4.5) and (4.1), one can obtain 


(4.4) 


(4.5) 

(4.6) 

(4.7) 


[(1 + £ 2 , r)(A 4 + E 7 )] t + — (1 + £ 2 t)(A 4 + A 7 ) < C£ 2 (A 4 + E 7 ) + C'5e 2 (l + £ 2 t) 

+ Cg J \0 y \(b\ + b 3 )dy + Ce(1 +£ 2 t)z j J -j^-\G\ 2 d£dy 

< C5£ 2 ( 1 + £ 2 t ) -1 + CK% + C(1 + £ 2 t)^K 5 . 


(4.8) 


Integrating (4.8) with respect to r and using (4.2) and (4.7), one has 


(A 4 + E 7 ) < CVS( 1 + £ 2 t)- 1+CoV *, / (1 + £ 2 s)(A 4 + K 7 )ds < CVS(1 + £ 2 t) C ov/J . (4.9) 

Jo 


Therefore, it holds 

||(0,V5C)(r)|| 2 < c^£ 2 (l + £ 2 t)- 1+Co ^. 

Multiplying (3.92) by (1 + £ 2 r)i _1? with fJ > (1 in Theorem 3.1 and using (3.14), (4.1), (4.7) 


(4.10) 
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(4.9) and the smallness of 6 , one has 


[(1 + e 2 T)l-*E 7 ]r = (5 - 0)(1 + £ 2 t)^£ 2 E 7 + (1 + £ 2 t) 3 2 -^E 7t 


< 


csk 8 + c(i + £ 2 t)I-»(k 4 + k 7 ) + C 5 (i + e 2 t) Y - 2 \\d^,C )\\ 2 

M=i £ 

+ C(1 + s 2 t)I-»Y 11^.^011“+ §(l + e 2 r)i- tf ||(^ > ^C)ll E H 9< W,0l| 3 

|a|=l |or|=l 

+ Ce (1 + e 2 r) 5 _1?+C ' 0%/ ^ J J ^Ml\G\ 2 d^dy + C 8 e 2 {1 + eV) -1- * 

< C 5 K 8 + C 5 { 1 + e 2 t)(IQ + K 7 ) + C(1 + £ 2 t)^ +CoV *K 5 + Cfe 2 ( 1 + eV) -1- * ( 4 . 11 ) 

+ ^ ( i +£ 2 T)f^ ^ h^, 0,011“ + §(1+6^)1^11(0,^011 E ii^(0,0,oii 3 - 

<4-i H=i 


By using (3.11) and (4.9), one can get 


\ j (i+ £ 2 s )2 * e 11^(0,0,0113^ + 4 f (i+ £ 2 s) 3 ^ i[( 0 ,0,on E iid Q ( < 0'OOii 3 d'S 

6 1/0 H = 1 6 1/0 |a|=l 

1 


/ (l + e 2 s)f- ,? [(l + £ 2 s)-i + i Co ^ + (l + e 2 s)- 1+Cov/J ] E l|d“( 0 , 0 ,Oll 2 ^ 


|a|=l 


1 


e J 0 


(l + e 2 s)5 ||9“(0,0,OH 2 ^ < Ce [ (l+£ 2 s)hQds <CVd£(l + £ 2 T) CoV *, (4.12) 


M=i 


provided that CoV^ < d. Thus integrating (4.11) over [0,r] and using (4.2), (4.7), (4.9) and 
(4.12) yield that 


E 7 < cVS{ 1 + £ 2 t)-^ +C o ^, 


which immediately implies 


(ik0^,oj 2 +Eii^^oii 2 +E / / 

l«l=2 M=2 J ^ 

+ - E [ [ ^W~ d ^ d y - C ^ 1 + e 2 r)-i +,?+Co ^. (4.13) 


|a|=l 


Proof of Theorem 3.1 Combining (4.3), (4.6), (4.13) and (4.10) and using the Sobolev 
inequality, it holds that 

!!(<*>, w) Iliac < C||(4>, T, W)|| (ll(0, 0, oil + fe(l + e 2 r)-2) < CVSe, (4.14) 


and 


\G[ 

M* 


/^4<^(l + £ 2 rri (4.15) 


Therefore, (4.6), (4.10), (4.13), (4.14) and ( 4.15| ) verify the a priori assumption (3.10) if we 
choose Ao = <5s. Hence, the proof of Theorem 3.1 is completed. □ 
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Proof of Theorem 2.4 The proof of (2.44) can be obtained directly from (3.9) by using the 
transformation (3.1) of the scaled variables (y. r) and the original variables {x, t). By combining 
(2.44) and Sobolev inequality, (2.45) can be derived immediately. Thus the proof of Theorem 
2.4 is completed. □ 
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